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QUASI-CYLINDRICAL THEORY OF WING-BODY INTERFERENCE AT SUPERSONIC SPEEDS AND 

COMPARISON WITH EXPERIMENT 1 

By Jack N. Nielsen 


SUMMARY 

A theoretical method is presented for calculating the flow 
field about a wing-body combination employing bodies deviat- 
ing only slightly in shape from a circular cylinder. If the 
combination possesses a horizontal plane of symmetry, no 
restrictions are required on wing plan form in the application 
of the method to the zero angle-of-attack condition. If the 
combination is lifting, the method requires that the wing lead- 
ing edges be supersonic. Then the extent of the flow field 
that can be calculated depends on the wing aspect ratio and 
whether or not the trailing edges are supersonic. Two methods 
of calculating the flow field, the W-function method and the 
multipole method, are presented. The methods as presented 
are accurate to the order of quasi-cylindrical theory. 

The method is applied to the calculation of the pressure field 
acting between a circular cylindrical body and a rectangular 
wing. These calculations are for combinations for which the 
effective aspect ratio of the wing panels joined together is greater 
than 2 and for which the effective chord-radius ratio is j or less. 
Two cases are calculated, the case in which the body remains 
at zero angle of attack while the wing incidence is varied and 
the case in which the wing remains at zero angle of incidence 
while the body angle of attack is varied. It was found that 
four Fourier components of the interference field are required 
to establish the pressure field, but that only one component is 
necessary to establish the span loading. A detailed discussion 
of the physical nature of the interference pressure field is given. 

An experiment was performed especially for the purpose 
of checking the calculative examples. The investigation was 
performed at Mach numbers of !Jf8 and 2.00 with a rectangular 
wing and body combination. Both the variable wing-incidence 
and angle-of-attack cases were covered. It was found that 
for sufficiently small angles, about 2° or less, the present method 
predicts the pressure distributions within About ±10 percent 
for both cases. Important nonlinear effects were found for 
angles of attack and incidence of 4° to 6°, and important 
viscous effects were usually found where laminar boundary 
layers encountered shock waves. 

INTRODUCTION 

In recent years the problems of supersonic wing-body 
interference have occupied the attention of many workers in 
aerodynamics. ■ The large amount of effort expended on the 
subject is a result of the important effects that interference 
can have on the overall aerodynamic characteristics of 


wing-body combinations. The trend toward using large 
bodies and small wings at supersonic speeds, especially for 
missiles, is the prime reason for the increased importance 
of wing-body interference at these speeds. 

Much significant work has already been done in the field. 
In reference 1, Spreiter has shown that when a wing-body 
combination is slender in the sense of his paper simple expres- 
sions for the lift and moment coefficients can be derived. 
These results were obtained by reducing a three-dimen- 
sional problem for the wave equation to a two-dimensional 
problem for Laplace’s equation. Another approach is that 
of simplifying the differential equation by using conical 
boundaries. Following this approach, Browne, Friedman, 
and Hodes in reference 2 obtained a solution for the pressure 
field of a wing-body combination composed of a flat tri- 
angular wing and a cone both with a common apex. The 
use of all-conical boundaries reduces the problem to one of 
conical flow for which powerful methods of solution are avail- 
able. 

Several investigators have presented methods for deter- 
mining the pressure field, including the effect of interference, 
acting on wing-body combinations employing circular 
fuselages and wings not necessarily slender. In reference 3, 
Ferrari has given an approximate method of obtaining the 
“interference of the wing on the streamlined body, assum- 
ing that the induced field generated by the wing is that 
which would exist around the wing if it were placed in the 
uniform stream alone.” Similarly, the interference of the 
body on the wing has been determined. The results of 
Ferrari thus represent a first approximation, and while a 
second approximation using the method is possible in prin- 
ciple it appears that too much labor would be involved. 
Morikawa in reference 4 has obtained an approximate solu- 
tion by solving a boundary-value problem and has also 
obtained a closed solution by approximating the three-di- 
mensional model by a planar model. Bolton-Shaw in refer- 
ence 5 has obtained a solution by satisfying boundary con- 
ditions at a finite number of points rather than over a surface. 

Another method for estimating the effect of interference 
on the aerodynamic properties of wing-body combinations 
which are not necessarily slender is given in reference 6. 
In this reference the method is applied to determining the 
drag of symmetrical wing-body combinations; it is also ap- 
plicable to the calculations of the lifting pressures acting on 
combinations employing wings with supersonic edges. In 


i Supersede? NAOA TN 2677 “Wing-Body Interference at Supersonio Speeds With an Application to Combinations With Rectangular Win©,” by Jadr.N, NWsen and William O. 
Pltt^ 1052, and NAOA TN 3128 "Comparison Between Theory and Experiment for Interference Pressure Field Between Wing and Body at Supersonio Speeds," by William 0. Pitts. Jaofc 
N. Nielsen, and Maurice P. GlonWddo, 19M. ./...• ,'r ■. > - • ! ’ 17 , 
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reference 7, an essentially new method of solving a wide 
class of wing-body interference problems has been presented. 
The method is based on decomposing the interference of a 
wing-body combination into a number of Fourier components 
and solving the problem for each component in a manner 
similar to that used by von Kdrm&n and Moore in reference 
8 for bodies of revolution. 

Phinney, reference 9, has compared the methods of refer- 
ences 3, 6, and 7 by applying each to the calculation of the 
pressure field acting on a circular cylinder intersected by an 
oblique shock wave. In reference 10 the theory of reference 
7 has been applied to the computation of the pressure dis- 
tributions acting on a rectangular wing and body combina- 
tion with the body at zero angle of attack and the wing at 
incidence. In reference 11 Bailey and Phinney have ap- 
plied the method of reference 7 to the calculation of the 
pressures on the body of a rectangular wing and body com- 
bination at angle of attack but with the wing at zero angle 
of attack. In reference 12 the same authors have com- 
pared their calculations with some experimental measure- 
ments made at a Mach number of 1.9. In reference 13 the 
experimental pressure distributions acting on a rectangular 
wing and body combination at Mach number 1.48 and 2.00 
are extensively compared with theoretical calculations based 
on the method of reference 7. 

In part I of the present report the theory of wing-body 
interference for combinations employing quasi-cylindrical 
bodies is presented, including recent developments not pre- 
viously reported in references 7, 10, or 13. The theory is 
applicable to combinations at zero angle of attack with 
horizontal planes of symmetry or combinations at angle of 
attack if the wing leading edges are supersonic. In part 
II the theory is applied to the calculation of the pressures 
and span loadings for a rectangular wing and body combi- 
nation for the case of the body at zero angle of attack and 
variable wing incidence and for the case of the wing at zero 
wing incidence and variable body angle of attack. The cal- 
culations for the second case are more complete than hitherto. 
In part III extensive comparison is made between the cal- 
culations of part H and the result of experiments at Mach 
numbers of 1.48 and 2.00 especially designed to check the 
calculations. 

SYMBOLS 
body radius, in. 

aspect ratio of wing formed by joining exposed 

half-wings together 

chord of rectangular wing, in. 

• c 

effective chord-radius ratio, — 

Pa 

strength of multipole of order 2 n at point x of 
body axis 

chord at wing-body juncture, in. 
chord at wing tip, in. 

arbitrary functions of s 

velocity amplitude function of nth Fourier com- 
ponent, in./sec 

wing-incidence angle, radians except where other- 
wise designated, positive for trailing edge down 
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modified Bessel functions of 
kinds, respectively 
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the first and second 


lift of combination back to wing trailing edge, 
lb; Laplace transform operator 
inverse Laplace transform operator 
lift on exposed half-wings joined together, lb 
lift on exposed half-wings in combination with 
body, lb 

index identifying sets of multipole solutions 
free-stream Mach number 
characteristic functions for obtaining multipole 
strengths 

number of Fourier component 
static pressure, lb/sq in. 
static pressure in free stream, lb/sq in. 
static pressure at any particular orifice of wing- 
body combination when a B =i w — 0, lb/sq in. 
static pressure at wind-tunnel wall orifice, lb/sq 
in. 


fQ rQ 2U 

pressure coefficient, - — —> — ^ for theoretical 
(Zo V 

calculations 

interference pressure coefficient due to nth 
Fourier component 

free-stream dynamic pressure, lb/sq in. 
dynamic pressure based on condition at wall 
orifice of wind tunnel, lb/sq in. 
cylindrical coordinates: y—r cos 6, z=r sin 0 
(See fig. 1.) 

Reynolds number based on wing-chord length 
real part 

semispan of wind-body combination, in.; Laplace 
transform of x coordinate 
axial, lateral, and vertical perturbation velocities, 
respectively, in./sec 
free-stream velocity, in./sec 
characteristic functions for calculating pressure 
coefficient 

Cartesian coordinates: x, axial coordinate; y, 
lateral coordinate; z, vertical coordinate, in. 
(See fig. 1.) 

body angle of attack, radians excopt whcro other- 
wise designated 

upwash angle of body-alone flow, radians 
wing angle of attack, radians 

VM*=I 

effective aspect ratio 
Dirac delta function; 

d(x)=0,x^0; f 3(x)cfo=l 

— co 

polar angle (See fig. 1.) 

taper, — 
c r 

dummy variable of integration 
sweep angle of wing leading edge 
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<P interference perturbation velocity potential 

<P i n nth Fourier component perturbation velocity 

potential 

(fa combination perturbation velocity potential 

<Pw wing-alone perturbation velocity potential 

(pi y R wing-alone perturbation velocity potential due to 

the exposed right half of the wing 
ipw L -wing-alone perturbation velocity potential due to 

the exposed left half of the wing 
tpw D wing-alone perturbation velocity potential due to 

the portion of wing inside the region occupied 
by body 

$ Laplace transform of <p 

SUBSCRIPTS 

L lower surface of combination 

U upper surface of combination 

I. GENERAL INTERFERENCE THEORY 

PHYSICAL PRINCIPLES 

Prior to a mathematical formulation of the wing-body 
interference problem, it is well to define interference and to 
explain how it arises. With a stationary wing or a stationary 
body in a uniform parallel flow, there are associated the 
wing-alone and body-alone flow fields. The wing-alone flow 
field does not, in general, produce flow tangential to the 
position to be occupied by the body surface. As a result 
an interference flow field must arise to cancel the flow field 
induced normal to the body by the wing. For this reason, 
the sum of the body-alone plus wing-alone flow fields will 
not be the flow field for the body and wing together. The 
difference between the flow field of the body and wing to- 
gether and the sum of the body-alone and wing-alone flow 
fields is defined to bo the interference flow field. 

The effects of wing-body interference on the flow field of 
a wing-body combination are illustrated by considering 
separately the effects of each component on the others. For 
the purposes of this discussion figure 1 shows a wing-body 
combination divided into the part in front of the leading 
edge of the wing-body juncture, henceforth called the nose, 
the winged part and the part behind the wing trailing edge, 
henceforth called the afterbody. If the combination pos- 
sesses a horizontal plane of symmetry and the angle of attack 
is zero, no restrictions on wing plan form are necessary. 
However, if the wing is twisted or cambered or if the nose 
is at angle of attack, then the wing leading edges must be 
supersonic for the following discussion to apply. 



Figure 1 . — Components of typical wing-body combination. 


Effect of nose on wing. — Consider now the flow as it 
progresses past the body. At the body nose the flow is that 
around a body of revolution, and it can be treated by existing 
methods such as those of references 8 and 14. When the 
body is at angle of attack a Bl there is an upwash field in the 
horizontal plane of symmetry of the body. If the body is 
sufficiently slender, the flow field in a plane at right angles 
to the body axis corresponds to that around a circular 
cylinder in a uniform stream of velocity, V sin a B . This 
gives an upwash field in the horizontal plane of symmetry 
of the body of 

a v =a B (l-\-a i /lf) (1) 

The effect of this upwash on the wing can be obtained by 
considering the wing to be at angle of attack and twisted 
according to equation (1) and by applying the formulas of 
supersonic wing theory. The wing pressure field so obtained 
is exact, within the limitations of the theory, for that section 
of the wing outboard of the Mach line emanating from the 
leading edge of the wing-body juncture. If the wing is 
located close to the body nose so that there is a chordwise 
variation in the upwash field due to the body, then the wing 
is effectively cambered, and the solution is more difficult. 
However, for most wing-body combinations it is possible to 
disregard the effect of the nose, and to assume that the wing 
is attached to a circular cylinder that extends upstream 
indefinitely. 

Mutual effects between body and wing. — The mutual 
interference between the body and wing on the winged part 
of a combination causes an interference field acting. on the 
body and on the wing inboard of the Mach line emanating 
from the leading edge of the wing-body juncture. The 
wing-alone flow field does not, in general, produce flow 
tangential to the position to bo occupied by the body sur- 
face. An interference flow field must arise that cancels the 
velocity induced by the wing-alone flow field normal to the 
body while not changing the wing shape. Alternately, the 
origin of the interference field can be explained in the follow- 
ing manner. The wing and body can be thought of as 
sources of pressure disturbances that radiate in all direc- 
tions in downstream Mach cones. The wing disturbances 
which radiate toward the body are, in part, reflected back 
by the body onto the wing and in part transmitted onto the 
body giving rise to interference pressures. Likewiso, the 
disturbances originating on the body pass onto the wing and 
affect the pressures there. It is apparent that the determi- 
nation of the interference pressure field on the body and on 
the wing inboard of the Mach line of the juncture is the crux 
of the wing-body interference problem. 

Mutual effects between wing panels. — To determine the 
region of influence of one wing panel on another, it is neces- 
sary to trace the path of a pulse from one wing panel across 
the body onto the other. The path traced across the body 
by the pulse originating at the leading edge of the wing-body 
juncture is the forward boundary of the region of influence of 
one wing panel on the body. (See fig. 1.) It is clearly the 
helix intersecting all parallel elements of the cylinder at the 
Mach angle. The boundary crosses the top of the body a 

distance of — 1 downstream and reaches the opposite 

wing-body juncture a distance downstream. A 
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pulse originating at a point on one wing panel and traveling 
to a point on the other panel can travel around the body on 
its surface to the opposite juncture and then along the wing 
to a given point, or it can leave the body tangentially before 
reaching the opposite wing juncture in a straight path to the 
point. The second means of transmitting the impulse is 
shorter in distance than the first and is the one which 
determines the forward boundary of the region of influence 
of one wing panel on the other. Applying this consideration 
to the pulse originating at the leading edge of one wing-body 
juncture, it is easy to show that the forward boundary of the 
region of influence of one wing panel on the opposite wing 
panel is given by the equation 



This boundary is also shown in figure 1, and it becomes 
parallel to the Mach line at distances far from the body. 

Effects on the afterbody. — As far as the interference 
effect of the body on the wing is concerned, it is confined to 
the winged part of the combination, but the effect of the wing 
on the body is felt also on the afterbody. For a symmetrical 
configuration at zero angle of attack there is no downwash in 
the horizontal plane of symmetry and the afterbody presents 
no particular problem. However, behind a lifting wing there 
is a downwash field. If the downwash were known every- 
where in the wing wake, then the wake could be considered as 
an extension of the wing with twist and camber. The wing 
wake and afterbody could then be incorporated with the 
winged part of the combination and treated in the same 
manner. However, the actual downwash pattern in the 
wing wake depends on the interference effect of the body on 
the wing. It is thus apparent that the solution of the after- 
body problem requires that the interference problem for the 
winged part of the combination be solved first. Only the 
winged part of the combination is analyzed in detail in this 
report. 

Regions of applicability of tbe theory. — The present 
interference theory can be applied to all or part of a wing- 
body combination depending on the configuration and the 
lift. If the combination is not lifting and possesses a hori- 
zontal plane of symmetry, then the interference pressure 
field can be determined for the entire combination. For a 
lifting combination with subsonic leading edges the upwash 
field in front of the wing makes the present method in- 
applicable. 

For a lifting combination with supersonic leading edges 
several geometric factors considerably influence the difficulty 
of calculating the interference field or indeed the extent to 
which it can be calculated. The effect of one of'these factors, 
the sweep of the trailing edge, is illustrated in figure 2 (a) . 
A subsonic trailing edge gives rise to multiple Mach wave 
reflections which greatly complicate the determination of 
the interference field over the rear part of the wing. Another 
important effect limiting the applicability of the theory is 
illustrated in figure 2 (b). This figure indicates that the 
interference field behind the incident wave can influence the 
tip upwash field which, in turn, influences the pressure field 
behind the reflected wave from the tip in a complicated way. 




(a) Subsonic trailing edge. (b) Effeot of -wing-body Intorfor- 

(c) Simple easo. ence on tip upwash Sold. 

(d) Tractable case. 

Figure 2.- — Classes of interference problems for lifting wing and body 

combinations. 

Avoidance of this complication requires that the incident wave 
intersect the trailing edge rather than the wing tip. To 
assure this condition, the aspect ratio must be greater than 
a certain Tninirnmn value in accordance with the following 
inequality: 

fiA> 7 -Q— - 77 - 

One of the simple cases of wing-body interference for a 
lifting wing and body combination is shown in figure 2 (c) . 
Here the leading and trailing edges are both supersonic, and 
the root-chord Mach wave intersects the trailing edgo. 
Also the wing-tip Mach wave intersects the body downstream, 
of the wing-body juncture so that no wing-tip effects occur 
on the wing interference pressure field. This condition 
imposes the aspect-ratio inequality: 

0A> A-: TV 

(i+x)(i+!Sp) 

Under the circumstances of this figure, the interference 
problem proceeds as if the combination had a horizontal 
plane of symmetry. Any body upwash field in front of 
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the wing can be treated as equivalent to a change in thick- 
ness distribution. The rectangular wing of aspect ratio 
greater than two is an example of the simple case, and it 
will ho treated as an illustrative example in this paper. 

An example of a tractable although fairly complicated 
case to which the present theory can be directed is shown in 
figure 2 (d). In region A the pressure field is determined 
as a pure wing-alone problem with any body upwash being 
treated as equivalent to a change in thickness distribution. 
In region B the problem is still a wing-alone problem which 
is complicated by upwash outboard of the tip. In region 
O there are body interference effects but no tip effects. 
In region D both effects prevail. In region E the tip has 
influenced the flow at the body surface and produced a 
secondary effect on the interference pressure field. 

MATHEMATICAL FORMULATION OF PROBLEM 

Throughout the analysis, the body radius is taken as unity 
and M? is taken as 2 so that >5=1. Any formula can be 
generalized to any body radius by dividing all length symbols 
by a, and to any Mach number by dividing all streamwise 
lengths by /3, by multiplying all pressure and lift coefficients 
by /3, and leaving all potentials, lift forces, and span loading 
unaltered. It is necessary to specify the wing alone before 
any detailed interference calculation can be carried out. 
However, in the theoretical solution of the problem the 
wing-alone definition is arbitrary. The flow field about 
the combination does not depend on the definition of the 
wing alone. 

General, decomposition of boundary- value problem. — 
The general case of a combination at angle of attack with 
the wing at incidence as shown in figure 3 is considered. 
The mathematical details of the decomposition of this 
configuration into tractable configurations is carried out in 
detail in Appendix A following the suggestions in reference 
16. A simplified discussion of the decomposition is now 
presented. The complete combination can be decomposed 
into three component configurations as shown in figure 4 (a) 
in which the wing boundary conditions are to be applied 
in the 2=0 plane and the body conditions on the r=l 


y y 



Figure 3. — General combination under combined effects of angle of 
attack and wing incidence. 


cylinder. Component (1) is simply the body alone, which 
creates an upwash field a„ in that region to be occupied 
by the wing in accordance with equation (1). Components 
(2) and (3) are combinations with wings of the same plan 
form; but while component (2) has a wing at angle of attack 
i w , component (3) has a wing with angle of attack — a u . 
The significance of this particular method of decomposing 
the general wing-body problem is that component (1), the 
body alone, can be solved by known methods and com- 
ponents (2) and (3) with bodies at zero angle of attack 
can be solved by the methods of this report. In the wing- 
incidence case where a B =0, only configuration (2) remains. 
This configuration can be decomposed into a wing-alone 
problem and a distorted-body problem as shown in figure 
4 (b). We co nfin e our attention to this wing-incidence 
case for the time being. 



$c $w + * 

(b) 

(a) Decomposition of general wing-body combination. 

(b) Decomposition for wing-incidence case. 

Figure 4. — Decomposition of wing-body combinations into simpler 

combinations. 

Consider now a combination with the body at zero angle 
of attack and let <pc be its potential. (See fig. 4(b).) This 
potential can be considered the sum of a wing-alone po- 
tential and of an interference potential <p. 

<Pc =< P TV~\~V (3) 
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Since the body is an infinite circular cylinder at zero. angle of 
attack, it produces no flow field. If the body were quasi- 
cylindrical with small distortions, a potential due to the 
body could be included in equation (3) . If the body has a 
horizontal plane of symmetry, the inclusion of a potential 
due to body distortion will not change the interference 
potential. 

The essential problem is to determine <p. First, select a 
convenient way of extending the wing through the body to 
form the wing alone, thereby specifying The wing-alone 

flow field in general produces velocities normal to the 

surface that will enclose the circular cylinder as illustrated 
in figure 4(b) for the region above the wing. In figure 4(b) 
and subsequent figures, all bodies are shown as cylinders 
parallel to the x axis. While the bodies of the component 
configurations in some cases are slightly distorted cylinders, 
they are nevertheless shown as true cylinders. This pro- 
cedure is compatible with the fact that the boundary con- 
ditions are to be applied on a true cylinder. The value of 

varies with 0 and with x. This means that a body con- 
forming to the wing-alone flow field is distorted in a compli- 
cated fashion. Now since the body must be circular, there 
must arise an interference potential <p that identically cancels 

^ 2 - at the body surface, thereby straightening it. 


c )<p cJ^ik 

dr dr 


at r=l 


( 4 ) 


There are two other conditions to be fulfilled by <p. It must 
not distort the shape of the wing when added to to pro- 
duce <pc . Thus when 0 = 0 , 


1^=0 
r d 0 U 


( 5 ) 


or iir =0 for the interference combination as shown in figure 
4(b). The last condition is that the interference potential 
must be zero ahead of the winged part of the combination. 


¥ 5 = 0 , x<0 


( 6 ) 


Equations (4), (5), and ( 6 ) are the essential boundary con- 
ditions on <p. 

The normal velocity ~ to be induced at the body surface 

by the interference potential can be analyzed at any given 
streamwise position as a Fourier cosine series. The ampli- 
tudes of the various Fourier cosine terms, / 2 „ (x), vary with x, 
the streamwise distance. Thus, 


^=^/ 2 *(a:) cos 2 n 0 =-^ 


dr t= 


at r=l 


( 7 ) 


Only even multiples of 0 are considered because of the vertical 
plane of symmetry. Consider that the interference po- 
tential is decomposed into a series of potentials such that 
each cancels one Fourier component of the velocity at the 
body surface; that is, 

<p2n (8) 

n— o 


with 

cos 2nd at r=l (9) 

Then the combination giving the interference potential $ can 
be decomposed in a series of combinations, each giving one of 
the <p 2 » values. The decomposition is illustrated in figure 6 . 



Figube 5. — Decomposition of interference combination into serios of 
Fourier component interference combinations. 


For n — 0 



and there is no variation of the normal velocity, pressure, or 
potential with 0 . Thus the first interference combination is 
a body of revolution. The pressure field acting on the body 
of such a combination can be determined by the method of 
reference 16. This n ~ 0 interference combination has the 
very simple significance that its flow normal to the r— 1 cyl- 
inder, subtracted from reduces the flow across the 
’dr dr 

body to zero when averaged from 0=0 to 0 =v at any stream- 
wise location. For 7 i=l, 

cos 20 

and the normal velocity, pressure, and potential will vary as 
cos 20 . 

To summarize briefly, it has been shown that the general 
interference problem of a body and wing at different angles 
of attack can be broken down into wing-body problems with 
bodies at zero angle of attack as shown in figure 4 (a). 
Combinations with the body at zero angle of attack are 
decomposed into wings alone plus interference combinations 
as in figure 4 (b). The interference combinations are finally 
decomposed into their Fourier components as in figure 6 . 

A general method for determining the characteristics of 
any Fourier component will now be given. It will be shown 
that good accuracy can be obtained for the interference 
potential with few Fourier components. 

SOLUTION BY METHOD OF IK FUNCTIONS 

The problem to be solved is that of a supersonic wing and 
body combination subject to the conditions already men- 
tioned, but with the wing and body possibly at different 
angles of attack. This problem is reduced to a body-alone 
problem and two wing-body problems with the body at zero 
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angle of attack as shown in figure 4 (a). The body-alone 
problem can be solved by existing methods such as references 
8 and 14. The procedure necessary to solve either wing-body 
problem as given in reference 7 is now summarized together 
with recent improvements. 

The potentials ipc, <pw, and ip must all fulfill the equation of 
linearized compressible flow 


( M a — l ) v ?«— < p n — <®«=0 (10) 

If we restrict ourselves for the time being to the case M= -yj2 
and transform equation (10) to polar coordinates, we have 

¥>rr+- 0 (11) 

with the coordinate system of figure 1. In solving the prob- 
lem we change from the physical space, x, r, 8, to the trans- 
formed space, s, r, 8, by means of the Laplace transformation 


independent of the boundary conditions. The inverse 
transform of the product of the two transforms can then be 
determined by the convolution integral. The part of the 
transform independent of the boundary condition can be 
thought of as defining a set of characteristic functions or 
influence coefficients. A tabulation of these functions allows 
a numerical solution of the problem for all boundary con- 
ditions. 

The manner of splitting equation (17) into two transforms 
depends also on the existence of the inverse transforms of 
the parts into which it is split. Let us write equation (17) ns 


cos 


2 ndFi (s)e~ ,(r ~ v . n ^ aa ^ i 1 1 

L &,'(*) 1 Jr Vd 


With the aid of the following relationships 

L~ l [F 3n (8)e-^]=Mx-r+l) 


(18) 

(19) 


/■» a> 

L [<p(x)\ e~ ai <p(x)dx=^(s) (12) 

With the boundary condition given by equation (6) that <p is 
zero 2 for x<0, equation (11) can be transformed to 

■fVr+i (13) 


L~ 1 (s^)=<p x (20) 

and the definition of the characteristic functions 

(2i) 

we obtain 


Expanding $ in a cosine serios of multiples of 8, we can 
satisfy the boundary conditions given by equation (5), and 
sinco there is a vortical plane of symmetry, we can co nfin e 
ourselves to even multiples of 8. With this restriction, gen- 
eral solutions to equation (13) can be written 


cos 2nd [ G ln (s)K la (sr)+D 2x (8)I 2n (8r )] (14) 

n=o 


where I 2x (sr) and K iH {sr) are modified Bessel functions. 
The constants C 3k (s) and Z) 2 ,(s) are arbitrary functions of s. 
The functions / 2rt (sr) can logically be eliminated at this point 
since from their asymptotic forms they can be shown to 
represent waves traveling upstream. The function C in (s) 
can bo evaluated by means of the remaining boundary con- 
dition given by equation (7). If we let 


then 


F 2H (s)=L[f ix (x)] 


c 3 m= 


sK 2n '(s ) 


(15) 

(16) 


Wo then have as the operational solution to our problem 


*=g cos 2n6F 2n (8) (17) 

The solution can be split into the product of two trans- 
forms, one dependent on the particular boundary conditions 
as represented by the F 3n (s) functions and another part 


i The condition frequently stated In deriving equation (13), that &,*=«0 for s—0 1 -. Is not 
required as proven in reference (17). This Is in accord with the Intuitive physical idea that 
any step in <£■ at the origin can be replaced by a continuous curve which for engineering 

purposes can have an effect different from that of the step only in a limited local region. 


<p t =± cos 2n0 F f* m~r+l)W 2a (x-^~ fUx i' +1) 
»■=■<> -yjr 

cos 2nd [ f r " r+ / 2 ^)W Sn (r— r+1— ^,r)d^- f2 " (a: ~ r+1) 

n«o | o T 

(22) 

With the aid of the W 2n (x,r) fimetions, the value of <p x , and 
hence the pressure or potential anywhere, can he calculated 
from equation (22) by numerical integration for as many 
harmonics as desired. This result was previously given (refs. 
7 and 10) for the r=l case only as 

cos 2 nd [" f j 3 n{Z)W in (x—Qdi—j 3K {x)\ (23) 

n-o L J o _I 

and the W 3x ( x) functions were tabulated for numerical cal- 
culation of the body pressure distributions only. The gen- 
eralization of the W 3n (x) function to IFa„(a;, r ) functions by 
means of equation (21) is a natural extension that permits 
the simple calculation of the pressure anywhere in the flow 
field. Some mathematical properties of W 2x (x, r ) functions 
and methods for their evaluation by automatic computing 
machinery have been studied by Dr. W. Mersman of the 
NACA. A r6sum6 of his results is reported in reference 18. 

Properties of the W 2n (x,r) functions. — Two important 
properties of the W 2n (x, r ) functions that make them useful 
for numerical work are that they possess no singularities in 
the field and their magnitudes never become large. These 
advantages are in distinct contrast to several disadvantages 
of the multipole method subsequently to he described. 
Curves of the W 3n {x, r) function are presented in chart 1 for 
»=0, 1, 2, and 3 for use in numerical computations. 
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A simple physical picture of the W 2n (x, r ) functions can he 
obtained from equation (22). Write the interference pres- 
sure coefficient 3 due to any harmonic as 


'2\ b<p 2 n _2 cos 2«$ r/a.Cc— r+1) 

v) dx~ v l ^ 


X x-r+l 

f2n(ZWsn(z—r + 1 



Let the velocity amplitude function be a delta function at 
the origin as shown in figure 6. Then 


_2cos2n « , r) " 


(25) 



MI 

v 


^llnit areo 

/ 


X 


Florae 6. — Fourier component interference combination with “delta 
function” protuberence at x=0. 


It is seen that physically the W 3h (x, r) function represents 
the pressure field due to a delta function in the velocity 
amplitude function. The first term represents an infinite 
pulse propagated along the Mach cone with apex at x= — 1 
and attenuating inversely as -Jr. The W iu (x— r+1, r) term 
represents the overexpansion behind the bump where the 
pressure would be zero if the flow were two-dimensional. 

Formulas for the W in (x, r) functions for small x and large 
x can be obtained from Laplace transform theory. In fact 
these results are 

(J*) (26) 


W SK (x,r)- 


4(r**-|-r 2 “)(47i+l)! 
' 2 in+1 (2nl) s 




(27) 


* The equation tor pressure coefficient to be used with equation (23) can on theoretical 
grounds be changed In going from the wing to the body of the configuration as discussed in 
reference 11. However, for simplicity, the linearized form of the B emouHi equation Is retained 
throughout the theoretical calculations. The contribution of the quadratic terms to the 
body pressure coefficient is subsequently discussed for the case of the combination at angle 
of attack. 


W( ^”i^( 3 +?)-^( 11+ 7+?)+°M » 

[(16n 2 +3)- (16 ^~ 1) ] (29) 

It should be noted that the value of W 2 „(0, r) is known 
precisely. 

SOLUTION BY METHOD OF MULTIPOLES 

Multipole types. — In references 7 and 10 a multipole 
method was used to determine the pressure field off the body. 
The singularities arising in this method, together with the 
loss of accuracy for the higher harmonics due to large numbers, 
led to the development of the W lK (x, r) method just de- 
scribed. Since the multipole method has application to 
certain problems and since its connection to tho W 3n (x, r) 
method is of interest, it will be given here. In the W ln ( x, r) 
method the pressure field'is determined by using boundary 
conditions on the body surface and continuing the pressure 
field outward from the body. It is intuitively obvious that 
any quasi-cylindrical flow can be generated by distributing 
sources and multipoles along the body axis in variable 
strength. If the strength of the axial multipole distributions 
can be related to the body shape (velocity amplitude func- 
tions), then the entire flow field can he calculated outward 
from the axis. 

Consider equation (14) which, -with D an (s ) equal to zero, is 


yift.fe - ) cos 2 nOK ta (sr) (30) 

n«=>0 

This equation can be interpreted to mean that the poten- 
tial is built up from a distribution of multipoles corresponding 
to the inverse transform of cos 2 ndK in (sr) along tho x axis 
in strength Ci n (x). However, there are many possible sets 
of multipoles corresponding to x integrals or derivatives of 
the set just mentioned. These are generated simply by 
rewriting equation (30) as 

[«"$„(«)] [ C092n ^» (gr) . ] (31) 

The first term represents the axial strength function, and the 
second term represents the fundamental multipole solutions. 
For each value of the index m there is obtained a distinct sot 
of multipole solutions. For selected values of m the multi- 
poles have the following forms: 


771 = 0 


COS 

771=1 
COS 


r, nT—irrs- t m cos 2n<? cosh (271 cosh-'r/r). ^ 
2ndL T-K 2 *(«■)]= p== — > x>r 

=0; x<y 

2 ML- [l^£i]_cos2«) x>r 

=0; a s<r 


( 32 ) 


771 = 277 
COS 


(z 3 — 7 •*)*"— 1/s ; 

=0; x<r 
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For positive values of m it is clear that at the Mach cone 
(x=r) no nonintegrable singularities occur. For negative 
values of m, derivatives of the m=0 multipoles are en- 
countered and the singularities occur on the Mach cone 
rather than on the axis. Since these singularities occur in 
tho flow field, they are not well suited to numerical methods 
of analysis. Another set of simple multipoles with singulari- 
ties on tho Mach cone are those of reference 19 given by 

r** cos 2nd 

( x S_ r 2)2»+ 1/S 


Multipole strengths. — The first step in determining the 
interference potential by the method of multipoles is to 
determine the multipole strength from the velocity amplitude 
functions. For the m= 0 set of multipoles the relationship 
between these two quantities is already given by equation 
(16) 

° 3 '^ = 8kL'1) (33) 

If the potential at a point P as shown in figure 7 is desired, 
the 


y,z 




Figure 7. — General point at which potential is to bo determined, 
multipoles must be distributed from —1 to x—r along the 
body axis. Since Laplace transforms must he zero for 
2 < 0 , the axial distribution must be shifted a distance at 
least unity to the right by introducing e~’ into the transform 

e ~ a0 ^= [^m\ [Fin{8)] m 

Equation (34) defines a Faltung integral involving a new set 
of characteristic functions given by 

M 2k {x)=L~ 1 (35) 

These characteristic functions have been studied in reference 
7, are tabulated in table I, and are plotted in figure 8. 
Forming now the 



Figure 8. — Graphical representation of functions. 


Faltung integral of equation (34) we obtain 

c,* (2-1)= JV*. (QM 2k (x-Q df (36) 

The function M 2x (x) has a square-root singularity at the 
origin so that Ci n (x— 1) will be finite if j 2n (%) is finite. How- 
ever, fin(x) may have a singularity which in confluence with 
the square-root singularity of M 2n (x) produces a singularity 
in Cinix— 1). 

Increasing the index m by unity has the effect of integrating 
the set of multipoles with respect to x and of differentiating 
the axial strength functions by x. While this decreases the 
order of the singularities of the multipole solutions, it in- 
creases the order of the singularities of the axial strength 
functions. The highest index m that does not lead to singu- 
larities thus depends on how many nonsingular derivatives 
Ca n ( x— 1) possesses, which, in turn, depends on the smoothness 
of fin(x)- In the calculations for the wing-incidence case 
(ref. 7) ftn(x) has a square-root singularity at 2=1 so that 
£ 2 /( 2 — 1) has a logarithmic singularity. Since <*,,,( 2 — 1) 
corresponds to m=0, it was possible to use multipole solu- 
tions of the m= 1 class and still obtain integrable singular- 
ities in the axial strength functions. 

Properties of the M 2n (x,r ) functions. — The M in (x,r) 
functions have simple physical significance. Let the velocity 
amplitude function he that corresponding to a delta function 
as shown in figure 6. Then by equation (36) 

c 2n (x—l)=M 2n (x) 


Thus the M 2n { x) function represents the distribution along 
the axis of multipole strength for the m — 0 set of multipoles 
necessary to make the velocity amplitude function a delta 
function. Correspondingly, it is the distribution necessary 
to produce apressure field corresponding to the W 2n {x-r-\-l,r) 
function. Equations (25) and (30) yield the relationship 
between the W 2n {x,r) and M 2n {x) functions. 


7 2»(x— r+1, r)=J^ M 2h (x—Q 


cosh ( 2n cosh 1 


Vft+l) 3 -ri 








Series for the M 2n { 2 ) functions for small and large values 
of the argument have been obtained by the standard methods 
of Laplace transform theory in reference 7 


102!£±> ' S+ 2W-3^+33 I 

7T l_VZ 4 JO J 

(38) 

The square-root singularity of M 2a (x) at the origin is note- 
worthy. For the asymptotic result only a single term has 
been calculated 

M 0 (x) 1 (39) 


MU x) 


— 16»(6a)! 
~ (27il) 8 2 8 * +1 



(40) 
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Span loading and pressure coefficient. — From equation 
(31) the potential can be written 


cos 27id[e~ , C 2n (s)] [K 3n (sr)e?] 


so that 


<P=jtj cos 2nd 




<hn(x— 1—?) cosh (2n cosh -1 


V 


r J 


V(£+i Y-t* 


(41) 


% (42) 


From this result the potential can readily be obtained and 
hence the span loading. The pressure coefficient follows 
directly from equation (42) using the linearized, form of 
Bernoulli’s equation 

P=- (43) 

Thus 

9 f C2„(— 1 ) cosh ( 2 n cosh -1 

COS Snfll „ , T -++ 

v n**o I 


y^+i)’-^ 

C* b Ca *( a:— ^ — ^) c°sh (in cosh -1 ^ 1 

Jr-l 55 V^+i) 2 -^ J 


(44) 


The practicability of using this result for calculating the 
pressures depends in the first place on the accuracy with 

which ~ C 2 *(x— 1—0 can be calculated. Since this calcula- 
tion depends on the M 3n (x) functions, which are tabulated 
at the present time only to the third decimal place, only 
three significant figures will usually he obtained for the axial 
strength functions. For higher harmonics and large values 
of x, loss of accuracy is incurred through the nature of the 
multipole solutions themselves. The following tabulation 
illustrates the point. 


cosh (2a cosh -1 x) 


X 

0 

1 

2 

3 

1 

1 

1 

1 

1 

2 

1 

7 

97 

1, 350 

3 

1 

17 

580 

19, 600 

4 

1 

31 

1920 

119, 000 


Although the set of multipoles used here is not well adapted 
to'the calculation of pressure coefficient for high harmonics 
and large values of x, it nevertheless is useful for calculating 
span loadings since only one or two harmonics are needed 
in this case. The difficulties of computing pressure coeffi- 
cient are alleviated in part by the fact that the pressure 
disturbances due to higher harmonics damp out within a few 
downstream radii. In reference 7 the pressure coefficients 
were computed up to the fourth harmonic (n=3) hut with 
some difficulty. The use of a set of multipole solutions other 
than the m=0 set does not hold much promise since increas- 
ing the value of m introduces singularities into the axial 
strength functions and decreasing the value of to introduces 
singularities into the multipole solutions. "While these sin- 


gularities are tractable using the methods of analysis, they 
are not adapted to the numerical methods used herein. The 
method utilizing the W in (x,r ) function avoids difficulties 
with large numbers and with singularities. 

IL APPLICATION OF THEORY TO COMBINATION OF 
CIRCULAR BODY AND RECTANGULAR WING 

In this part of the report calculations are carried out to 
determine the pressure field acting on a wing-body combina- 
tion employing a rectangular wing with no thickness. The 
calculations are first made for the body at zero anglo of attack 
with the wing at incidence — the wing-incidence case. The 
calculations are then made for the body at angle of attack 
with the wing at zero angle of incidence — the angle-of-attack 
case. For the calculations attention is focused on the upper 
half of the combination since the experimental measurements 
were made for the upper half. 

WING-INCIDENCE CASE 

The complete pressure field will now be calculated. As 
previously mentioned, the wing alone can be specified in 
any convenient manner and, for the purpose of the example, 
the wing alone is taken as the rectangular wing extending 
straight through the body from side to side. Although the 
analysis as carried out is for M=-yJ 2, the results are pre- 
sented in a form applicable to a range of Mach numbers. 
The steps in performing the calculation are: (1) to determine 
tpw, the wing-alone potential; (2) to determine the velocity 
amplitude functions, f 3n (%); and (3) to determine the 
potential or pressure, as desired, anywhere in the field. No 
tip effects are considered until the results are presented as a 
function of wing aspect ratio. 

Wing-alone potential. — The wing-alone flow, exclusive 
of tip effects, can be determined from the Ackeret theory. 
The flow at a spanwise station out of the region of influonco 
of the wing tips is illustrated in figure 9. The potential for 
the flow above the wing is 


<Pw= Vx when z>x (45) 

<Pw=Vx+i w V(x—z) when z< x (46) 

The sidewash produced by such a potential is zero, and 


/ 



Figube 9. — Flow field produced by wing alone. 
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, /x FiTrr2cos(27i-l) w 2 cos (271+1) w 

/ s .( x )=— 2^=1 


271+1 

4 


4 n- 


y 


when a;<l 


when a; - 1 


(51) 

(52) 


Figure 10. — Variation of normal velocity induced at body surface by 
wing alone; unit body radius; wing-incidence case. 



Figure 11. — Graphical representation of velocity amplitude functions; 
wing-incidence case. 


the downwash is uniformly —i w V. The down wash 
oz 

causes a flow normal to the surface r=l in amount —iwV 
sin 6. This means that, for a body conforming to the wing- 
along flow, the deformation is zero at the wing-body junc- 
tures and a maximum on the top of the body. The inter- 
ference combinations when added to the deformed body 
straighten it out Fourier component by Fourier component. 

Fourier amplitudes of body normal velocity. — The Fourier 
amplitudes of the normal velocity induced by the wing-alone 

potential at the body are determined by expanding at 


r=l in a Fourier cosine series of even multiples of 8. The 
normal velocity distribution is shown in figure 10. For 
£>1 the body is totally immersed in the wing downwash 
field. With the usual equation for obtaining the Fourier 
amplitudes of a function, there is obtained 


o peln-ix 

f 0 (x)=- i w V sin 8 dB (47) 

"ft J 0 

A (*skr l x 

fn s (x)=- iwV sin 8 cos 2nd dd (48) 

ttJo 


The integrations give 

» fo(z)=— — (l— VI— a?) when x<l (49) 

7T 

f 0 when x>l (50) 

t r 


where w=sin -1 x. The/ 2 „(a:) functions are shown in figure 11. 
The constant values of / 2 „( x) for »>1 are noteworthy. 
The values of / 2 *(a;) are tabulated in table II. 

Interference pressure distributions. — The interference 
pressure distributions have been calculated for the first 
four Fourier components and are presented in figure 12. 
In this figure the abscissa is proportional to distance behind 
the Mach line originating at the leading edge of the juncture, 
as illustrated in part (a) of the figure. Although the calcula- 
tions have been carried out for M—-^2, that is, /3=1, and 
for unit radius, they are generalized to all Mach numbers 


X 7* 

and body radii by replacing x— r+1 by ^— -+ 1 and P 2 „ by 

p P 2 » as has been done in the figure. From the figure It is 
apparent that the cusps in the pressure distributions are 

X T 

propagated downstream along lines of constant ^ — -+1 or 

x—pr; that is, along the downstream characteristics. As 
the pressure distributions move outward from the body 
along the downstream characteristics, they are distorted 
and decreased in magnitude. 

Increasing the order of the Fourier harmonics causes two 
important effects: first, the number of points of zero pressure 
is increased and, second, the pressure coefficient damps 
more rapidly. As a result of the first effect, the contribu- 
tions of the higher harmonics to the combination span 
loading are proportionally less than their contributions to 
the pressure coefficient; while, as a result of the second 
effect, the more remote a point is from the leading edge of 
the wing-body juncture, the fewer the number of Fourier 
components that must be included to obtain its pressure 
coefficient accurately. All interference pressure distributions 

x t‘ 

exhibit discontinuities in slope at -= 1-1=1- This be- 

pa a 

havior is a consequence of the fact that the body becomes 
totally imm ersed in the wing-alone flow field for this condi- 
tion. When the pressure distributions of the various Fourier 
components are added together to obtain the interference 
pressure distributions, the discontinuities in slope tend to 
cancel so that the pressure distribution for the combination 
will be smooth. 


A detailed examination of the interference pressure dis- 
tribution for the first Fourier component illustrates several 
points of interest. The importance of the component arises 
from the fact that it accounts for most of the effect of inter- 
ference on the span loading. The reason for this is that the 
pressure coefficients for n = 0 are of invariable sign. The 
effect of the first Fourier component is to reduce the velocity- 
induced normal to the body by the wing-alone flow field to 
zero averaged around the body from 8 = 0 to 0=ir at any stream- 
wise location. 


'■'fy & ez s ° 3 *i/ z d $ 
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0 .5 LO 1.5 2.0 2.5 30 3.5 4.0 


x/fio - r/o + I 

(a) n= 0 

(b) n=l 


Figube 12. — Interference pressure distributions of various Fourier components; wing-incidence case. 
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5 2D 

x/fia-r/o +1 

(c) n=2 

(d) n=3 

Figure 12. — Concluded. 
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For purposes of comparison with the exact results for 
7i=0, some approximate results have been included in 

figure 12 (a). For values of j^~0 on the body, the Ackeret 

value of P 0 (twice the local stream angle divided by /9) is a 
close approximation to the true pressure coefficient. This 
is the result of the facts that the part of the body affecting 
the interference is effectively plane for points near the lead- 
ing edge of the wing-body juncture and that there is no 
variation of any quantities with 6 so that an approximate 

two-dimensional situation prevails. As — increases beyond 

pa 

unity on the body, there is a rapid decrease in the pressure 
coefficient below the Ackeret value due to the effect of all 
disturbances in front of the point in question as represented 
by the integral of equation (24). 

In reference 7, the following approximate results were 

obtained for small and large values of ■£- — -+1 for the pres- 
to- a. 

sure coefficient: 


pPiK ~^(&ri +1 ) cos 2n0 


(53) 


X r l t n 

as 3 J— 1 — >0 

/3a a 




aD . 32 i w (4ra!)(r 2 “+r- J «) cos 2nd / x \~ «»+« 

P ln 7r(2a!)W— 1)2 4 “ +1 \j3a/ 


(54) 


x 

as — — > oo 
pa 

% T 

For -r hi <0.6 equation (53) is a good approximation for 

pd d 

n=0 although it is of little value for higher-order harmonics. 
There is a general tendency of P 0 to approach a unif orm value 

independent of r as — 


r 

a 


-1 becomes large, as shown by 

equation (54). The damping in the characteristic direction, 
although initially inversely proportional to the square root 
of r, is ultimately independent of r. 

Pressure distribution in juncture of wing-body combina- 
tion. — By adding the interference pressure coefficients of 
the various Fourier components to that for the wing alone, 
the pressure distribution for the combination is obtained. 
The addition has been carried out for the wing-body juncture 
using four Fourier components and six Fourier components, 
and the results are presented in figure 13. The pressure 
coefficient with interference is less in magnitude than 2, the 
value without interference, showing that significant losses of 
lift occur in the wing-body juncture. A comparison of the 
results for four components and six components shows that 
four components give good over-all accuracy for all values of 


^greater than 1. 


For small values of — in the wing-body 


juncture, the curvature of the body insofar as the flow is 



Figure 13. — Theoretical pressure distribution at wing-body juncture 
of combination using four and six Fourier components; wing- 
incidence case. 


concerned is not large so that the body is effectively a vortical 
boundary on which a given distribution of normal velocity 
is producing an interference field. Supersonic wing theory 
applied to this condition gives for the net interference pres- 
sure coefficient (ref. 7). 

pP 4x , a; .... 

iw Sir Pa Pa v ' 

It is clear that the calculated results can be jointed smoothly 
to this result. Using the result of equation (56) onablos 
satisfactory results to ho obtained with four Fourier com- 
ponents. 

The critical region in the convergence of the solution is 
that near the leading edge of the wing-body juncture. The 
higher harmonics have their most important effect near hero 
and rapidly damp downstream along the body. Hence more 
and more Fourier components would be required to got 

*c 

accuracy for smaller and smaller values of However, 

pa 

with the result of equation (55), this extra work is unneces- 
sary. 

One point of interest in figure 13 is the fact that when — 

pd 

equals approximately 3, the pressure coefficient increases in 

X 

magnitude. This is due to the fact that for— >ir the influ- 

pa 

ence of the opposite half- wing is felt in the wing-body 
juncture. 

Pressure distribution on top meridian of wing-body 
combination. — The pressure distribution on the top meridian 
of the wing-body combination is obtained in the some fashion 
as that at the wing-body juncture, the difference being that 
the pressures due to the even number Fourier components 
have the same sign at the meridian as at the juncture, whereas 
the odd numbered components have reversed signs. The 
pressure distributions based on four and six Fourier com- 
ponents are shown in figure 14. 

Several interesting effects are exhibited by the results. 

• • x 

The step in the wing-alone pressure at — =1 is effectively 

pd 

canceled by the interference pressures of the Fourier com- 
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FiauRE 14. — Theoretical pressure distribution on top of combination 
using four and six Fourier components; wing-incidence case. 


XX X 

ponents from — =1 to — =*-/2 , and for — >7r/2 the pressure 

increases rapidly and tends toward the two-dimensional 
value. The effect of the interference pressure in canceling 
the effect of the wing alone on the top of the body from 

^=1 to —=7i-/2 is to be expected since the wing of the 

combination can have no effect on the top of the body unless 

fls has been already pointed out. If an infinite 

number of Fourier components had been taken, the pressure 


• XX 

coefficients would be identically zero from — =0 to 

The general behavior in this regard is evidence of the plausi- 
bility of the calculated results. 

The tendency of the pressures to approach an asymptotic 
value is also illustrated by figure 14. This asymptotic 
value represented by the sum of the wing-alone pressure plus 
the asymptotic results foi the first Fourier component is 
given by the following equation: 


PP 

iw 


o 

- — ~ — Z- 


irx/Pa 


(56) 


For — )>2.4, the results of this equation are in good 

agreement with the results of figure 14. 

Some evidence is furnished from the pressure calculations 
for the juncture and top of the body concerning the number 
of Fourier components necessary for accuracy. Comparisons 
made in figures 13 and 14 show that about four components 
are sufficient and that the addition of two more is not worth 
the extra work. 

Pressure distribution on wing of wing-body combination. — 
The distribution of the pressure acting on the wing of the 
combination can he determined in a manner similar to that 
for the wing-body juncture by adding to the wing-alone 
pressure those due to the Fourier components. The result- 
ant pressure distribution for the wing based on four Fourier 


• • X 

components is shown in figure 15. For small values of — 

pa 

the higher-order oscillations in the pressure coefficient as 
shown in figure 13 have been ignored, and the curves have 
been faired through them. 



x/fi a 

Figure 15. — Theoretical pressure distribution acting on wing of 
combination; wing-incidence case. 


Since the region of influence of the body on the wing is 
confined to the wing region downstream of the Mach lines 
emanating from the leading edge of the juncture, in front 
of this line the pressures are uniform at the two-dimensional 
value, and behind the line there is a decrease in the magnitude 
of the pressure coefficient. If the body were a perfect 
reflector, that is, a vertical Avail of infinite extent, then 
there would be no pressure loss. However, the pressure 
pulses originating on the wing are only in part reflected by 
the circular body. The efficiency of the body as a reflector 
is discussed subsequently in connection Avith span loading. 
The tendency of the pressure to increase in magnitude near 
the inboard tra il ing edge is due to the effect of the opposite 
wing panel which at the Aving juncture is felt doivnstream 

of the point — =7r. 

Pa 

Span loading. — The span load distributions for a range of 
rectangular Aving-body combinations Avith the body at zero 
angle of attack can be determined from the pressure distri- 
butions of figures 13, 14, and 15. Since the pressure dis- 
tributions of figure 15 are in a form independent of Mach 
number, it is convenient to define a span loading which is 
an integral of these distributions. 


ga w ar 


=/:.[ne -§>©>( id 
-ic&regM* 


(57) 


The quantity in the square brackets is taken to be the span 
loading. If all distances are taken in units of the body 
radius, then “a” can be set equal to unity in the formulas. 

The pressure results of figure 15 are for values of the effec- ’ 
tive chord-radius ratio of 4 or less and for values of the 
effective aspect ratio of 2 or greater. Span loadings for 


c 

any combination of — (or c*) and pA in these ranges can be 

obtained by integrating the pressure distributions. The 
span loading evaluations have been made for c*=4 and 
pA>2. First the span loadings due to the- various Fourier 
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components are discussed, and then the span loadings for 
the actual wing-body combinations are presented. 

In figure 16, the contributions to the span loading for the 
first three Fourier components are shown. For n = 0 the 
pressure field does not depend on 0, being axiallyfsymmetric, 
and a constant loading exists on the body. However, on 
the wing as the spanwise distance increases there is a decrease 
in the span loading, due primarily to decrease in the length 
of chord over which the interference pressures act. The 
span loading due to the first Fourier component causes a 
loss of lift everywhere along the span. 



Fioube 16. — Theoretical span loading of various Fourier components 
acting on combination of body and rectangular wing having effective 
chord-radius ratio of 4; wing-incidence case. 


A comparison of the results of figure 16 for n=0 and 7i=l 
shows that the first Fourier component accounts almost 
entirely for the effect of interference on the span loa ding of 
the combination. For the body this fact is even more true 
than for the wing. This fact is of considerable importance 
since it gives a simple means of extending the lift and moment 

£ 

results to larger values of — than those for which the pressure 

pd 

distributions have been calculated. Also, it suggests a simple 
means of minimising the adverse effects of interference on 
lift as will subsequently be pointed out. 


With the techniques of Laplace transform theory, it is 
possible to obtain asymptotic formulas for the span loadings 
of the various Fourier components. For the first Fourier 
component the following asymptotic result has been obtained 
by the standard methods of Laplace transform theory. 
(See Appendix B.) 


- 2 f (f ) * (SHrffi * - 1 ** (3m)+ik 


when — — > to 
P a 


(58) 


The asymptotic result for the span loading given by this 
equation, when compared with the results of the exact cal- 
culations in figure 1 6, is seen to be slightly low. However, for 

values of — greater than 4 the difference between the results 

pd 

decreases, and equation (58) thus provides a satisfactory 
means of extrapolating the results of the present calculations 

for span loading to larger values of — 

pd 

The asymptotic result has also been determined for the 
higher-order Fourier components as a matter of interest. 
The span loading is 




8 cos 2nd 


im(4n 2 — l)r* B 

16 (471—1) ! cos 2nd 


w(2rr!) a (47r a -l)2^- 1 (^ 


as - — >oo 
Pa 


( 60 ) 


The results of equation (59) and the exact solution for 7i= 1 
in figure 16 both corroborate the fact that the span loadings 

of all but the first Fourier component are negligible for ■§- >4. 

pd 

It is also to be noted that the contribution to the loading of 
the first component given by equation (58) increases without 
limit as x—> <» ; whereas the span loadings of the higher-order 
components are finite. 

To obtain the span loading for the family of combinations 
for which — = 4, it is necessary to consider the loadings of 

pd 

both the wing alone and the Fourier components. The 
necessary calculations have been carried out, and the span 
loadings for the family of combinations based on one and 
four Fourier components are both shown in figure 17. The 
loading due to the wing alone is also shown. No effect of 
wing tips has been included. It is to be noted in figure 17 
that, whereas the loading on the wing due to its own pressure 
field is constant, there is some loss on the body because of 
the fact that the pressure field of the wing alone acts on the 
body only if x>pa sin d. However, if an afterbody is 
included, some of the lift lost can be recovered. As has 
already been pointed out, the pressures due to the first 
Fourier components are positive on the upper half of tlio 
wing-body combination and produce a loss of lift, as figure 17 
shows. When the effects of four Fourier components are 
taken into account, the net lift is slightly higher than that 
for one Fourier component, but the difference is not signifi- 
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Fiotob 17. — Theoretical span loading for combination of body and 
rectangular wing having effective chord-radius ratio of 4; wing- 
incidence case. 

cant. For most engineering purposes, one Fourier compo- 
nent is sufficient for determining the span loading when 



Some insight into the mechanism of wing-body interference 
can be gained by comparing the span loading for the com- 
bination with those for two reference loadings: (1) the com- 
plete reflection case for which the blanketed area of the wing 
acts effectively at i w , and (2) the no -reflection case for which 
the blanketed area of the wing is supposed to act effectively 
at zero angle of attack. The span loading corresponding to 
the first case of complete reflection of the wing pressure pulses 
by the body is, in fact, the span loading marked “wing alone” 
in figure 17. A comparison of this curve with that based on 
one or four Fourier components shows that the loading given 
on the assumption that the wing blanketed area is fully effec- 
tive in lift is too optimistic. Under the conditions of the 
second reference loading, the sole purpose of the blanketed 
area is to support lift generated by the wings. A comparison 
of the span loading for this case with the true loading shows 
that the average load on the body is well predicted, but that 
the loading on the wing is underestimated. - A comparison of 
the true loading with those for the two reference cases reveals 


the interesting fact that the body is somewhat less than 50 
percent effective in reflection for this particular family of 
configurations. 

lift. — For values of — < 4 the pressure distributions 

fJw 

already presented are sufficient for obtaining span loading or 
lift on either the wing or body for all combinations having 
sufficiently large aspect ratios to avoid effects of the tips on 
the wing-body interference. This is the case for f3A>2. The 
lift results are presented in terms of a nondimensional param- 
eter k w , defined as the ratio of the lift on the exposed half- 
wings in combination (exclusive of that on the body) to that 
on the exposed half- wings joined together. 



O £ B =0 


(60) 


For — >4 the value of k w can be obtained by using the 
asymptotic form of the span loading given by equation (58) . 



c*+c* log (^pj)-J-log (c*+l) 
vd*(?pA— 1 ) 


c 


CO 


(61) 


The values of k w have been determined from the pressure 

£ 

distributions of figure 15 for values of — <4 and from equa- 

pd 

£ 

tion (61) for values 4 of — >4. The effect of the wing tips 

has been taken into account by utilizing reference 20. The 
results are shown in figure 18 wherein k w is given as a function 



Figure 18. — Lift effectiveness for wing or control surface in 
combination with body. 


of — for various effective aspect ratios of 2 and greater. It 
pa, 

should be borne in mind that the results of the figure are for a 
combination of body and rectangular wing or an all-movable, 
' rectangular control surface with no gap. It is noted in the 

c • • 

figure that the exact results for — <4 can be faired into the 

pa, 


4 In reference 10 the asymptotic analytical expressions for Jtirand I. VC are not precise by 
virtue of an incorrect upper limi t on an integral. The maximum numerical error in kw Is 
about 0.01 and in x^je about 0.03 for very large values of e*. The precise values are given 
in this report. Those for ffA—2 are unchanged. 
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• c • • 

asymptotic results for — >4, thereby providing a design 

pd 

chart for engineering purposes for the entire range of — 

fjd 

• c 

The curves of figure 18 illustrate the decrease of k w as — 

fXL 

increases at constant effective aspect ratio, and the slow in- 
crease of Jc w as the wing chord becomes very large. The loss 
of lift is most serious for pA=2, being about 15 percent in 
the worst case. 

A practical point in connection with the loss of lift on the 
wing due to interference is that this loss occurs no matter 
what the body angle of attack, even though the calculations 
are made for a B =0. It occurs either in the case of a wing 
mounted on a body or in the case of a deflected all-movable 
control surface. For wings with swept leading edges for 
which all of the wing area lies in the region affected by the 
interference, even larger losses than occur with rectangular 
wings are to he anticipated. However, the loss of lift at the 
design condition can, at least in principle, he largely pre- 
vented by designing the fuselage so that it conforms to the 
first Fourier component in the wing-alone flow. This would 
involve contracting the fuselage above the horizontal plane 
of symmetry in a rotationally symmetric fashion and expand- 
ing a like amount beneath the horizontal plane of symmetry. 
Whether or not such a change would improve the lift-drag 
ratio can best he determined by experiment. 

Center of pressure. — The center-of-pressure locations 
have been calculated for the same range as the lift results of 
figure 18. The center-of-pressure location in chord lengths 
behind the leading edge are presented in figure 1 9 . For large 
values of c*, an asymptotic result has been calculated for 
x C p/c using the methods of Laplace transform theory and con- 
sidering only one Fourier component. 


11 2 


— low 

( 2c* \ 

— | log (c*+l)] 

x ep 2 3 PA PAc * 3 

L *• + 

2 0g 


c 

k w 

fi 1 > 


0 2 pAj 


as c*-> <» 


(62) 



Figure 19. — Center of pressure for wing or control surfacp in 
combination with body. ■ 


The values of have been determined from the pressure 


distributions of figure 15 for values of -^-<4 and by equation 

pa 


(62) for values of — >4. 

pa 


The loss of lift near the tips has 


been taken into consideration. The exact results for — <4 

Pa~ 

have been faired into the asymptotic results for large values 
of — by dashed curves to provide an engineering design 


chart covering the entire range of — ■ 


It is again mentioned 


that this chart is applicable both to the wing of an airplane 
or missile or to an all-movable, rectangular control surface 
with no gap. The curves of figure 19 start at values of 

c 


X c 

corresponding to those for the wing alone at — : 


=■0. As 


pa 


increases for constant p A, there is a forward movement of the 
center of pressure because of the loss of lift due to interference 
which is mostly effective on the rear of the wing. For the 
lowest effective aspect ratio of 2 there is about a 4-percont 
forward movement of the center of pressure due to interfer- 
ence in the extreme case. For largo effective aspect ratios 
the forward movement is not nearly so largo. As the value 

of — increases for constant pA, there is an asymptotic 
pa 

approach of the center of pressure back to the wing-alone 
value. 

ANGLE-OF- ATTACK CASE 


In figure 4 (a) it is shown how the flow field of a combina- 
tion can he built up of a body alone and two wing-body flow 
fields. The first wing-body flow field ((2) of fig. 4 (a)) has 
been solved in the preceding section, and we now solve the 
second wing-body problem ((3) of fig. 4 (a)). The •wing is 
effectively twisted so that the slope of its surface is a u as given 
by equation (1) . It should be noted that the problem of the 
combination and body ■with a rectangular wing twisted 
according to the second term of equation (1) has been 
solved by Bailey and Phinney in reference 11 using the 
present theory. Their calculation is restricted to the body 

It 



Figure 20. — -Graphical representation of velocity amplitude functions; 
■ • . . angle-of-attack case. . 
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(a) n= 0 

Figure 21. — Interference pressure distributions of various Fourier components; angle-of-attack case. 


and is carried out for downstream distances of 2pa from the 
wing leading edge. Actually, the results of reference 11 
represent the difference between the angle-of-attack and 
wing-incidence cases treated here and are in agreement with 
the present results. This agreement is, in effect, an inde- 
pendent check on the accuracy of the present numerical 
results for the interference pressure distributions. 

Wing-alone potential. — The first step in the calculation 
is to determine the wing-alone potential. Because the wing 
is twisted to conform to the body upwash field, this determi- 
nation is fairly tedious and has been carried out in Appendix 
O. The form of the wing-alone potential found in reference 
11 is in agreement with those found herein for the wing- 
incidence and angle-of-attack cases. 

Fourier amplitudes of body normal velocity. — The velocity 
amplitude functions for the present case were computed 
numerically by performing a Fourier analysis of the calcu- 
lated body normal velocity distribution at a number of body 
cross sections. An analytical determination was made of the 
velocity amplitude functions by the authors of reference 11 
for values of x<2. However, for ri>2 the velocity amplitude 
functions are said by these authors to lead to incomplete 

413872—57 88 


elliptic integrals, and no analytical determination was made. 
A numerical determination has been made herein for 0<jr<(4:. 
The numerical values of the/ 2s (a:) functions for this case axe 
tabulated in table II and plotted in figure 20 for illustrative 
purposes. 

Interference pressure distributions. — The interference 
pressure distributions have been calculated by numerical 
integration using equation (22). The results are shown in 
figure 21. The interference pressure distributions are very 
similar to those for the wing-incidence case, being about 
twice as large. 

Pressure distribution in juncture of wing-body combina- 
tion, — The pressure distribution of the combination is 
obtained by adding the interference pressure coefficients to 
the pressure coefficients of the wing alone. The results, 
using four and six P 2 „ components, are shown in figure 22. 
This figure shows that four components give a close approxi- 
mation to the linear-theory value for xft 3a>l. At x/pa=0, 
the wing leading edge, linear theory with Beskin upwash 
theory gives exactly jSP/a= — 4.0. For the region r//9a<(l 
the higher harmonics have their greatest importance, and 
many components would be necessary to get good accuracy. 
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Figure 21. — Concluded. 



Figure 22. — Theoretical pressure distribution at wing-body juncture 
of combination using four and six Fourier components; angle-of- 
attack case. 


However, satisfactory accuracy can be obtained by fairing a 
curve through this region since both end points are known. 

One item of interest in figure 22 is the increase in the 
magnitude of f) P/a B near point 1. This is due to the influence 


of the opposite half-wing reaching the -wing-body juncture 
at this point as shown in the sketch. 

It should now be noted that the pressure distribution for 
the-ease of the body at angle of attack with the iving at zero 
angle of incidence is represented by the sum of cases (1) and (3) 
as given by figure 4 (a). However, we are neglecting tho 
contribution of case (1) because it is small. The contribu- 
tion to the pressure coefficient represented by case (1) is 
that due to a yawed infinite cylinder since we are neglecting 
nose effects. This contribution, which is clearly present in 
front of the wing, is 

(J)=a B (l-4cos ! 0) (63) 

For the juncture of the combination (0=0°) the contribution 
is about 0.1 for a B = 2° and 0.3 for a fl =6°. At a D = 2° tho 
effect is thus negligible compared to P/a s of about 4, and at 
00=6° there are definite nonlinear effects that make a pre- 
cise application of linear theory inaccurate. For those 
reasons the contribution given by equation (63) has boon 
neglected. For the top and bottom of the body tho contri- 
butions are one-third of the foregoing and hence are also 
negligible. 
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Pressure distribution on body of wing-body combination. — 
The pressure distribution on the body is also obtained by 
adding the interference pressure coefficients to the pressure 
coefficients due to the wing alone. The interference pressure 
distribution for any value of 9 differs from that in the wing- 
body juncture, 6=0, only by a cos 2 n6 factor. For example, 
in the juncture cos 2nd is always +1. On top of the body, 
0=7r/2, cos 2nd alternates between +1 and — 1 as n increases. 
On the 0=x/4 meridian cos 2nd has values of 0, +1, and —1, 
so that when n is odd P 2 ,=0. The pressure distributions on 
the top meridian of the body and on the 5=45° meridian of 
the body are shown in figures 23 and 24, respectively. 



Figure 23. — Theoretical pressure distribution on top of combination 
using four and six Fourier components; angle-of-attack case. 



x/pa 


Figure 24. — Theoretical pressure distribution on 9=45° meridian of 
body of combination using four and six Fourier components; angle- 
of-attaok case. 


Several interesting effects are exhibited by figures 23 and 
24. The step in the wing-alone pressure at x/fia=l in figure 
23 is effectively canceled by the interference pressure from 
x/f3a= 1 to x/Pa=ir/2, and for a;//3a>7r/2 the pressure increases 
rapidly. The effect of the interference pressure in canceling 
the effect of the wing alone on the top of the body from 
ar//9a=l to a;//3a=7r/2 is to be expected since the wing of the 
combination can have no effect on the body in front of the 
Mach helix (point 1 of sketch) originating at the leading edge 
of the wing-body juncture. If an infinite number of com- 
ponents had been computed, the combination pressure 


coefficients would be identically zero from x/pa=0 to x/pa= 
tt/2. The same effects are exhibited by figure 24 except that 
the wing-alone step occurs at x/fla=-fi/2 and the Mach 
helix intersects the meridian at a://}a=ir/4, point 1. The 
Mach helix from the opposite wing panel intersects the 
meridian at point 2 causing an additional pressure rise. 
Since the region in which fiP/a B =0 is known and since the 
exact linear theory is well approximated by four components 
for large values of x/fia, theoretical curves of good accuracy 
can be faired from figures 23 and 24. The area under the 
high peaks in the curves near »//3a==ir/4 would become in- 
finitesimal if an infinite number of interference pressure 
components were taken. 

Pressure distribution on wing of wing-body combination. — 
For the region in front of the Mach wave from the leading 
edge of the juncture, the calculation of pressure coefficients 
is just a wing-alone problem. The pressure coefficients 
in this region can therefore be obtained directly from the 
wing-alone potential as given in Appendix C. The result is 

^ >= 1 (® 4 ) 

In the region behind the Mach wave the pressure coeffi- 
cients were obtained directly from the W 2n (x, r) functions, 
as was done on the body. The results of these calculations 
for the wing pressure distributions are shown in figure 25 
and are to be compared with the pressure distributions of 
figure 15. 

Span loading. — The span loadings have been determined 
by graphical integration of the pressure-distribution curves 
of figures 21 to 25. For a combination with a value of 
c/j 3a of 4 the span loadings associated with the various 
Fourier components are shown in figure 26, which is to be 
compared with figure 16 for the wing-incidence case. The 
magnitudes for the n=0 harmonic of the angle-of-attack 
case are about twice those for the wing-incidence case, but 
otherwise the two cases are similar. The span loading 
including wing-alone and interference effects is shown in 
figure 27, which is to he compared with figure 17. The 
important difference is noted that the peak span loading is 
nearly equal to the root loading in the angle-of-attack case, 
but is considerably greater than the root loading in the wing- 



Figbre 25. — Theoretical pressure distribution acting on wing of 
combination; angle-of-attack case. 
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y/a 


Figure 20. — Theoretical span loading of various Fourier components 
acting on combination of body and rectangular wing having effective 
chord-radius ratio of 4; angle-of-attack case. 

incidence case. Because of the difference in the shape of 
the span loadings, a different trailing vortex pattern 'would 
be associated with each. No effect of wing tips is included 
in figures 26 and 27. 

lift. — From the theoretical wing pressure distributions of 
the combination the lift of the wing panels in the presence 
of the body can be calculated as a function of fiA and c/fia. 
To show how the body upwash is effective in increasing the 
lift of the wing, a factor K w has been calculated. This 
factor has been defined as 

iir = 0 (65) 

A TV 

Hero Live is the lift of the panels in the presence of the body 
and L w is the lift of the wing panels joined together at angle 
of attack a B . In calculating Lwc first the lift of the exposed 
panels as part of the wing alone must be calculated. This 
was done by the use of reversibility theorems described in 



Figure 27. — Theoretical span loading for combination of body and 
rectangular wing having effective chord-radius ratio of 4; anglo-of- 
attaok case. 


reference 21. The lift of the entire wing alone including 
the blanketed area was so determined. The lift of the 
blanketed area was then calculated from the potential 
function given in Appendix C and subtracted from the lift 
of the entire wing alone to get the lift of the panels. The 
loss of lift on the panels due to interference as determined 
by graphical integration was then subtracted to get L^ vc . 
The values of K w so calculated are shown in figure 28 (a) 
as a function of c/fia and in figure 28 (b) as a function of a/s. 
Figure 28 (a) shows a large effect of / 8 A at constant c/fia; 
whereas figure 28 (b) shows a small effect of j 8 A at constant 
a/s. 

In figure 28 (b) the effect of aspect ratio on K w at a fixed 
value of a/s is less than the precision of the calculations as 
indicated by the cross-hatched area. For comparison the 
values of K w calculated from slender-body theory have been 
included in the figure 



( 66 ) 
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Figure 28. — Lift effectiveness for rectangular wing in combination 
with body; angle-of-attack case. 



(b) Effect of radius-semispan ratio. 
Figure 28. — Concluded. 


The dose agreement between the linear-theory results for 
the present case and the slender-body-theoiy results is 
noteworthy since the rectangular wing and body combina- 
tions considered here are not slender. This result suggests 
that slender-body theory can be used for calculating lift 
ratios for nonslender configurations. 

m. COMPARISON OF EXPERIMENT AND THEORY FOR 
RECTANGULAR WING AND BODY COMBINATION 

APPARATUS AND PROCEDURE 

An investigation to evaluate the present theory was made 
in the Ames 1- by 3 -foot supersonic wind tunnel. This wind 
tunnel was equipped with a flexible-plate nozzle that could 
be adjusted to give test-section Mach numbers from 1.2 to 
2.2. The pressure measurements are obtained as photo- 
graphic recordings of a multiple-tube manometer board 
using dibutyl phthalate as the fluid. 



The sting-supported model, which is diagrammed in figure 
29, is a combination consisting of a cylindrical body with an 
ogival nose and a rectangular, wedge-shaped wing. The 
dimensions of the model are given in figure 29. The wing 
was made 10 percent thick to minimize aeroelastic effects. 
It was mounted in the body by means of a set of angle blocks 
which enabled the flat wing surface containing the orifices 
to be set at 0°, —1.9°, —3.8°, and —5.7° angles of incidence 
with respect to the body center line. The pressure orifices 
were all located on the upper surface of the model. The 47 
orifices were distributed along seven spanwise stations in 
order to give a comparison with theory for the wing and the 
body. The locations of the orifices are given in table m. 
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Since this investigation required a comparison of the data 
for several Mach numbers and Reynolds numbers at the same 
values of a B and i w , it was necessary to set a B and i w ac- 
curately for each measurement. The sting support by which 
the model was mounted had sufficient flexibility that it was 
deemed necessary to have means for accurately setting the 
values of a B and i w under tunnel operating conditions. The 
values of i w were accurately set by means of angle blocks in 
the body. The angle of attack was set by a special image 
projection device. A mirror was inserted in the scblieren 
system so that an image of the model was cast upon a screen. 
With the wind off, the model was set at the desired value of 
a B and the inclination of the model image was marked on the 
screen. With the tunnel in operation at the desired pressure) 
the angle of attack of the model was adjusted until the 
inclination of its image was parallel to the calibration line 
made on the screen with the wind off. To check this 
method, a horizontal and vertical wire grid was placed on the 
tunnel window and schlieren pictures were taken of the model 
while the tunnel was in operation. These pictures showed 
that the image projection device set a B to within ±0.07° of 
the desired value. It was especially necessary to set a B 
accurately for the small angles to avoid large percentage 
errors in the angle setting. 

The model angle of attack ranged from +6° to — 6° in 2° 
increments, and the wing-incidence angle ranged from 0° to 
—5.7° in 1.9° increments. The test was performed at the 
two Mach numbers 1.48 and 2.00 and at the Reynolds 
numbers of 0.6, 1.2, and 1.5 million, based on the wing chord. 
The model was tested for all combinations of these values of 
the four parameters investigated. 

A complete set of data in the form of P for the Reynolds 
numbers 0.6, 1.2, and 1.5X10® at Af=1.48 and for R= 
1.5X10® at Af=2.0Q,is presented in table IV. These values 
of P are, for the most part, averages of two readings. 

REDUCTION AND ACCURACY OF DATA 

All data are reduced to the coefficient form (p—pi)/g 0 - 
Actually the quantity (p— p r )/2r was measured, and subse- 
quent corrections were applied to change the reference static 
pressure to p v (pi is the static pressure at the particular orifice 
in question when a B =i w =0°) and the reference dynamic 
pressure to g 0 . Since p 1 includes the effects of nose thickness 
and stream angle, using p x as a reference pressure minimizes 
these effects and essentially gives only the pressures due to 
the angle settings of the model. The dynamic pressure was 
adjusted from q T to q 0 on the basis of a previous pressure 
survey of the tunnel. This latter adjustment was negligible 
for M— 1.48 and amounted to less than a 3-percent correction 
for M=2.00. For the purpose of comparison with theory 
the pressure coefficient ($—pi)/g 0 is reduced to the parameters 
P P/a B for ‘£ I r=0° and pP/i w for a B = 0°. 

Two types of errors entered into the experimental investi- 
gation: systematic errors and random errors. In this paper 
accuracy will be taken as the ability of the experiment to 
give the true values without nose effect or stream angle and, 
hence, is a measure of the systematic errors. Precision will 
be taken as the ability to repeat the data and, hence, is a 
measure of the random errors in the experiment. 


Several factors contributed random errors. The major 
factor was the error in the angle-of-attack setting. The un- 
certainty in each angle setting was ± 0.07°, but each measure- 
ment was dependent upon two angle settings: the setting for 
the condition represented and the setting to determine the 
zero correction. This leads to a net uncertainty of 0.1° 
which would account for a 5-percent error for angles of ±2°. 
Most of the remainder of the uncertainty in the data is duo 
to the fact that the reference wall static pressure in the tunnel 
changed slightly from run to ran while the total pressure 
remained constant. Although the magnitude of this pressure 
change was quite small, it was large enough compared to the 
small pressure differences for the 2° angle settings to cause 
as much as a 3-percent error. In addition to these factors, 
between 1-percent and 2-percent uncertainty was observed 
in reading the data from the manometer-board pictures. 

To determine experimentally the precision of the data, a 
large number of repeat measurements were taken and 
compared. It was found that for a B or i tr = ±2°, two 
independent determinations of pP/a B or pP/i w differed from 
each other by ± 7 percent on the average. For a B or i w = ±4° 
and a B or i w = ± 6°, the experimentally determined precisions 
of pP/a B and pP/iw are ±4 percent and ±2 percent, respec- 
tively. The precision in pPia B increases with the magnitude 
of the angle because a large part of the random error is duo 
to the angle setting. The known major experimental errors 
are due to stream-angle and body-nose effects. The effect 
of these factors was not determined, but, as previously 
described, corrections were applied to minimize their effect, 
assuming the effects did not vary appreciably with angle-of- 
attack settings. This assumption should bo good for the 
body-nose effect. However, it is not necessarily a good 
assumption for the stream-angle effect since the stream 
angle varies with vertical location in the tunnel and the 
model moves approximately 6 inches in a vertical direction 
between a B =4-6° and a B =— 6°. Since the stream-angle 
correction that was used was obtained for the a B =0° position 
in the tunnel, data obtained at a B =0° should have no 
appreciable error due to stream angle. For other values of 
a B , some error due to stream angle is possible.' 

For the purposes of this paper, the important question is, 
“How well does theory predict the experimental data?” 
Direct comparisons between linear theory and experiment 
will be made only for <^=±2° and i w =~ 1.9° data. In 
figure 30 experimental pressure distributions in the wing-body 
juncture obtained from two independent measurements with 
iw— — 1-9° and a s =0° are shown together with a faired 
curve of their average values. The ± 7-percent limit of 
precision about the average value is represented by the 
dotted lines. The figure shows that the theoretical value 
generally lies between these dotted lines, and therefore tho 
theory predicts the experimental values within the precision 
of the data in this example. 

GENERAL PHYSICAL PRINCIPLES 

Before the discussion of the results of the investigation in 
detail, it is well to give first a general physical description 
of the effects to be expected. Figures 31 and 32 show 

* A stream-angle and pressure survey of the wind ttmnel In the vertical plane of symmetry 
Indicated that stream-angle variation caused the magnltudo of tho experimental values of 
pPJaa to be 4 percent high on the average. 
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Figube 30. — Comparison between two independent readings of pressure 
distribution in wing-body juncture; a B = 0, *V=1.9°, Af=1.48, 
B«=ri.6X10<>. 


Inlersection of Mach 
cones with surface 
of combination 



Figure 31. — Isometric drawing of pressure distribution acting on 
combination of body and rectangular wing; angle-of-attack case. 

qualitatively the pressure distributions to be expected on a 
rectangular wing and body combination for the angle-of- 
attack case and the wing-incidence case, respectively. The 
cbordwise variations of the coefficient, fiP/a B or pP/i w , are 
shown for five stations by the shaded areas. 6 These figures 
show that Mach cones emanating from the wing-body 
juncture determine the points at which the various effects 
of wing-body interference are felt. On the cylindrical body 
the pressure coefficient is zero in front of the Mach helix 
originating at the leading edge of the wing -body juncture. 
The body pressure coefficients here are taken as zero because 
the effects of crossflow on the body pressures are very 
small, as shown in connection with equation (63) . However, 
as shown by the two stations on the body, the pressure rises 
abruptly behind this Mach helix, point 1, in both figures. 
The Mach helices from the two wing panels cross the 0=7r/ 2 
station simultaneously so that there is only one large increase 
in the magnitude of the pressure coefficient. These Mach 
helices cross the 0=3ir/4 station at two different points so 
that beyond point 1 there is a secondary increase in the 
pressure coefficients at point 2. These Mach helices continue 

« The pressure distribution shown for the 0=3x/4 station on the body Is Identical to the 
pressure distribution for the S<=xH station due to the symmetry of the model 
413072—57 84 


to curl around the body until they strike the wing panel at 
points 3, where part of the pressure disturbance continues 
along the wing and part of it is reflected along another 
Mach helix on the body, causing a further increase in the 
magnitude of the pressure coefficients at points 4. Another 
pressure disturbance originates at the trailing edge of the 
wing-body juncture that causes the decrease in the magnitude 
of the pressure coefficients noted at points 5 of the two 
figures. 


Intersection of Mach 

cones with surface 
of combination 



Figure 32. — Isometric drawing of pressure distribution aoting on 

combination of body and rectangular wing; wing-incidence case. 

On the wing of the combination the pressure coefficient is 
the same as that for a wing alone in front of the Mach wave 
from the wing-body juncture, except that when the body is 
at an angle of attack the body upwash effectively twists the 
wing in a manner such that a w =a B (l+a i /y !l ) . Figure 31 
shows this effect of body upwash along the leading edge of 
the wing where the pressure coefficient decreases as y/a 
increases because of the effective twist of the wing. The 
importance of body upwash can be seen by comparing the 
pressure distribution along the leading edge in figure 31 with 
that in figure 32. The pressure coefficient at the wing-body 
juncture in figure 31 is twice that in figure 32 where there is 
no body upwash. The pressure coefficient at any given 
spanwise station remains nearly constant between the wing 
leading edge and the Mach wave from the wing-body junc- 
ture. Behind the Mach wave, interference from the wing- 
body juncture causes the pressure coefficient to decrease in 
magnitude as shown in the two figures. 

EFFECTS OF ANGLE OF ATTACK 

Comparisons between theory and experiment for the 
angle-of-attack case are made in figures 33 for data at a 
Keynolds number of 1.5X10 6 and Mach numbers of 1.48 
and 2.00 with 7,^=0° and a B = ±2° and ±6°. 

Pressure distribution in juncture of wing-body combina- 
tion. — A comparison between linear theory and experiment 
for the pressure distribution in the wing-body juncture is 
made in figures 33 (a) and 33 (b) for both Mach numbers. 
The sketches show the pertinent Mach lines and the span- 
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(a) Af=1.48, y/a— 1.02 

(b) Af=2.00, y/a=1.02 

Figure 33. — Pressure distributions due to angle of attack; 12=1.5X10*. 

wise location of the orifices. 7 The experimental data points 
from the wing surface on which a compression occurs (nega- 
tive angle of attack) are represented by flagged symbols, and 
the data points from the surface on which an expansion 
occurs (positive angle of attack) are represented by un- 
flagged symbols. The figures show that the theory predicts 
the magnitude of pP/a B about 5 percent below the average 
of the a B = ±2° experimental values at M=1.48 and about 
15 percent below experimental values at M=2.00. The 
chordwiso variation is well predicted by the theory. 

Linear theory predicts that the parameter f} P/a B is inde- 
pendent of angle of attack. Actually it is not, and the 
nonlinear effects of angle of attack cause a spread in the 
data. It is possible to evaluate approximately the variation 
in the parameter fiP/a B with angle of attack at the wing 
leading edge. First the upwash just in front of the leading 
edge was calculated using equation (1) which is based on 
linear theory. Then the pressure coefficients at the wing 
leading edge were computed using shock-expansion theory. 
The values of pPja B for a B =— 6° and +6° so calculated are 

1 The location, of these Mach lines Is only qualitative because the calculations vrere made 
us Inc shock-expansion theory, with the assumption that there was no local Mach number 
variation behind the leading edge of the wing. To simplify the sketches, the Mach helices 
on the body are represented as straight lines. 


shown in figure 33 for values of y/a of 1.92, 2.58, and 3.92. 
For Af=1.48 body upwash caused the shock wave to bo 
detached from the wing in the wing-body juncturo so that 
no calculation of the spread could be made there. For 
M= 2.00 it was found that near the wing-body juncturo the 
predicted spread in |3 P/a B between —6° and +6° was about 
twice the experimental spread; whereas for y/a greater than 
about 1.5 the experimental spread was fairly well predicted. 
This difference between shock-expansion theory and tko 
experimental data in the wing-body juncture is probably duo 
to the combination of several things. First, near the wing- 
body juncture the body upwash is modified by viscous effects. 
Second, the theoretical spread was calculated at the leading 
edge of the wing, and this value was assumed to apply rear- 
ward to the first orifice. This assumption is probably good 
beyond y/a= 1.5 where the chordwiso changes in pressure 
are small back to the first orifice, hut, in the juncture, the 
changes in the chordwise direction are large near the wing 
leading edge so that this assumption is probably invalid. 
Third, the contribution of the body crossflow field previously 
mentioned is present (eq. (63)). 

Another phenomenon not predicted by linear theory is 
shown by figure 33 (a). The linear theory predicts that the 



(o) Af=1.48, top meridian, 
(d) Af=2.00, top meridian. 

Figure 33. — Continued. 




QUASI-CYUINDRICAL THEORY OP WING-BODY INTERFERENCE AT SUPERSONIC SPEEDS 


1327 



(e) Af=1.48, 0=45° meridian. 

(f) ilf=2.00, 0=45° meridian. 

Figtjrb 33. — Continued. 


Mach helix from the opposite wing panel (see sketch) should 
intersect the wing-body juncture at point 1, causing an 
increase in the magnitude of pP/a B . This effect is observed 
experimentally for negative values of a B in front of point 
1 rather than exactly at point 1. The reason is that for 
negative values of a B a compression occurs on the orificed 
surface reducing the local Mach number from the free- 
stream Mach number, thus increasing the Mach angle and 
causing the Mach helix to shift forward. The result is the 
spread of the data shown in figure 33 (a) near point 1. This 
effect is not shown by figure 33 (b) because the Mach helix 
lies more rearward for ilf=2.00 so that the orifices do not 
extend to the Mach helix as shown by the sketch. 

Figures 33 (a) and 33 (b) show that Mach number has no 
effect upon the magnitude of the higher-order spread due to 
angle of attack or upon the chordwise variation of 0 P/a B , 
but on the average the magnitude of 0 P/a B is about 10 
percent higher for M= 2.00 than for M=1.48. 

Pressure distribution on top meridian of body of wing-body 
combination. — A comparison between the linear theory and 
experiment for the pressure distribution on the top meridian 
of the body is made in figures 33 (c) and 33 (d). These 
figures show that theory and experiment are in good accord 



(g) M=1.48, j//o=1.25 

(h) M=2.00, y/a= 1.25 

Figure 33. — Continued. 


for a B = ±2°, particularly at Af=1.48. However, nonlinear 
effects due to a B cause a large spread between the data for 
a s =+6° and-a B =— 6°. All the effects predicted to occur 
on the body in the section of the report "General Physical 
Principles” are observed experimentally, but not exactly 
at the points predicted because of nonlinear effects. The 
pressure rise predicted at point 1 of figures 33 (c) and 33 (d) 
occurs prematurely and is less abrupt than expected for all 
angles of attack because of the boundary layer on the body. 
The variation in local Mach number causes the Mach 
helices to shift forward for the negative angles of attack 
as discussed in the section treating the wing-body juncture. 
The increase in the magnitude of 0 Pja B expected at point 2, 
x/0a=3v/2, actually occurs at about x/0a=4 for a B =— 2°. 
The decrease in magnitude of pP/a B that is expected at 
point 3 actually occurs at about x/0a-=4.O for a s =— 6°. 
For the positive angles of attack the Mach helices are shifted 
rearward so that these effects are not observed experimentally 
in the range of x/fia measured. 

Figures 33 (c) and 33 (d) show that, in general, the 
M=1A8 data are predicted better by the theory than are 
the M=2.00 data. For M= 2.00 there is an unexpectedly 
large pressure coefficient in front of point 1 for negative 
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x/Po 

(i) Af= 1.48, y/a= 1.92 

(j) Hf=2.00, y I a— 1.92 

Figure 33. — Continued. 

angles of attack. The predicted pressure coefficient due to, 
crossflow is only about 0.1 in units of the ordinate and hence 
does not account for the observed effect at M= 2.00. For 
a B = — 2° and M—2.00, pP/a B dips slightly near point 1 and 
then rises and overshoots the a B =— 6° data. This effect is 
due to the boundary-layer condition on the body and will 
be discussed in detail in the section de alin g with Reynolds 
number effect. 

Pressure distribution on 0=45° meridian of body of 
wing-body combination. — A comparison between the linear 
theory and experiment for the pressure distribution on the 
0=45° meridian of the body is made in figures 33 (e) and 
33 (f). Essentially the same effects are shown on this merid- 
ian as on the top meridian. 

Just as for the top meridian of the body the experiment is, 
in general, better predicted by the theory for Af=1.48 than 
for -M=2.00, and the same boundary-layer effects are evident 
near point 1 for M=2.00. 

Pressure distribution on wing of wing-body combination. — 
Experimental chordwise pressure distributions on the wing 
are shown in figures 33 (g) to 33 (n) for the four spanwise 
orifice stations y/a= 1.25, 1.92, 2.58, and 3.92. In front of 
the Mach cone from the wing-body juncture no interference 



x/$a 

(k) M=\M, yla= 2.58 

0) Hf=2.00, y)a= 2.68 

Figtjbb 33. — Continued. 

is felt from the wing-body juncture so that the theoretical 
pressure distribution for a wing alone in the body upwash field 
is used in this region. Figures 33 (g) to 33 (n) show that on 
the average the wing-alone theory predicts magnitudes of 
pP/a B about 5 percent below the measurements for a B — ±2° 
for .M=1.48 and about 12 percent below the measurements 
for M= 2.00. The spread in the data between a B =+ 6° 
and cc B — — 6° is fairly well predicted by shock-expansion 
theory for y/a greater than about 1.6 (figs. 33 (i) to 33 (n)). 
At yla= 1.25 the predicted spread (not shown) is too largo, 
just as for the wing-body juncture. 

Some of the interference effects discussed in the section of 
the report entitled "General Physical Principles” are illus- 
trated in figures 33 (g) to 33 (n). The interference effect from 
the opposite wing panel is observed in figure 33 (g) whero, 
just in front of point 1, the same spread in the data occuis as 
in the wing-body juncture. According to linear thooiy the 
disturbance originating at the nearer wing-body juncture 
should be felt at point 2 of figures 33 (i) to 33 (m), and the 
magnitude of ft P/a B should begin to decrease from the wing- 
alone value there. These figures show that the magnitude 
of /3 P/cc B does decrease in the neighborhood of point 2. They 
also show that, in general, the a n =-\- 6° and the a B —— 6° 
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(m) M=1.48, j//a=3.92 

(n) M=2.00, j//o<=3.92 

Figure 33. — Concluded. 


data come together in the neighborhood of point 2. This 
convergence is due to a variation in the local Mach number 
with ct D ■ This is shown by the sketch in figure 33 (j) where 
the disturbance from the wing-body juncture is first felt at 
point 3 for a B — — 6°, whereas it is first felt at point 4 for 
a n — +0°. Since the magnitude of f}P/a B begins to decrease 
as soon as this disturbance is felt, the magnitude of pP/a B 
begins to decrease at a smaller value of xjfia for a B = — 6° 
than for a fl =+6°, thus causing the convergence observed. 
The sketches in figures 33 (k) and 33 (m) show that the 
disturbance from the wing tip should also cause the a B = +6° 
and a D =— 6° data to come together beyond point 6 in these 
figures. The figures show that the data not only come 
together but actually cross over and reverse order just beyond 
point 6. 

Tho only significant effect of Mach number shown by 
figures 33 (g) to 33 (n) is the approximately 10-percent-larger 
values of pP/a B for M=2.00 than for M=1.48. Nearly 40 
percent of this difference may be due to differences in stream 
angle in the wind tunnel for the two Mach numbers. 

Span load distribution. — Span loa ding is defined for both 
the body and tho wing as the integral (see eq. (57)) 


The experimental and theoretical results for the span loading 
distribution on the wing and body of the combination are 
presented in figure 34. No account has been taken of tip 



(a) 2H=1.48 

(b) M= 2.00 

Figure 34. — Span load distributions due to angle of attack; 

R= 1.5X10’. 

effects in calculating the span loading because the twist of 
the wing makes a determination of these effects a difficult 
wing problem. The theory is thus valid only inboard of 
point 2. If an approximate answer is needed, the Busemann 
tip solution (ref. 20) can be joined onto the span loading at 
point 2. Figure 34 shows that the theory is generally about 
10 percent below experiment. This result is not surprising 
in view of the comparisons between the experimental and 
theoretical pressure distributions of figure 33. Of particular 
interest is the fact that, in general, the higher-order differ- 
ences due to a B that were so large for the pressure-distribution 
results are negligible for the span loading distribution. The 
only exception is on the top of the body, y/a=0, and M=2.00, 
where the effects of boundary-layer and shock-wave inter- 
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action are largo. The explanation for the independence from 
a B is that the higher-order effects on the top surface are 
compensated for by higher-order effects of the same magni- 
tude on the lower surface so that the net loading per unit 
angle is very nearly independent of angle of attack. 

EFFECT OF WING-INCIDENCE ANGLE 

Comparison is made between theory and experiment for 
the whig-incidence case for data taken at a Reynolds number 
of 1.5X10 6 and Mach numbers of 1.48 and 2.00 with a B =0° 
and i jr= — 1.9° and —5.7°. It will be remembered from the 
section on the accuracy of data that there is no appreciable 
orror duo to stream angle for the wing-incidence case, and 
the comparison between experiment and theory reflects this 
fact. 



(a) Ilf =1.48, y/a= 1.02 

(b) Af= 2.00, t//a=1.02 

Figure 35. — Pressure distributions due to wing incidence; 12=1. 5X10’ 

Pressure distribution in wing-body junctures. — The linear 
theory and experimental pressure distributions in the wing- 
body juncture are compared in figures 35 (a) and 35 (b). 
The symbols in the figures are flagged to be consistent with 
the use of flagged symbols for negative angle-of-attack data. 
The figures show that the experimental values are about 5 
percent below those predicted by the theory for i w =— 1.9°. 
The magnitude of the nonlinear effects due to i w is predicted 
at the leading edge by shock-expansion theory. Figures 



x/fla 


(c) Af=1.48, top meridian. 

(d) M =2.00, top meridian. 

Fioure 36. — Continued. 

35 (a) and 35 (b) show that the spread predicted in this 
manner can account for the experimental results. The pro- 
mature increase in the magnitude of pP/iw near point 1 is 
due to the effect of the opposite wing panel and variation of 
the local Mach number as discussed in the anglo-of-attack 
section. No significant effect of Mach number was found 
on the parameter pP/i w - 

Pressure distribution on top meridian of body of wing-body 
combination. — A comparison between the linear theory and 
experiment for the pressure distribution on the top meridian 
of the body is made in figures 35 (c) and 35 (d). Thoso 
figures show that theory and experiment are in good accord 
for i„r=— 1.9°. However, nonlinear effects duo to i w causo 
much larger differences between theory and experiment for 
i w = — 5.7°. This is consistent with the angle-of-attack 
case where the higher-order effect due to a B was large for 
negative angles of attack. 

All of the effects observed for the angle-of-attack caso duo 
to disturbances from the wing are also shown to occur for tho 
wing-incidence case in figures 35 (c) and 35 (d). Tho paths 
of these disturbances as predicted by linear theory are shown 
on the sketch, and the positions at which the effects are 
expected to occur are shown on the abscissa. 
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(o) A/=1.48, 0=46° meridian, 
(f) 1U=2.00, 0=45° meridian. 

Figure 35. — Continued. 


(g) Af=1.48, y/a= 1.25 
(K) M—2.00, y/o=1.25 

Figure 35. — Continued. 


The only significant effect of Mach number apparent in 
figures 35 (c) and 35 (d) is the larger boundary -layer and 
shock-wave interaction for M— 2.00 than for M= 1.48 near 
point 1. The M=2.00 experimental data for iw= — 1.9° 
dip and then overshoot at this point. This phenomenon is 
discussed in more detail in the section of the report on 
Reynolds number effect. 

Pressure distribution on 0=45° meridian of body of 
wing-body combination. — Linear theory is compared with 
experimental results for the pressure distribution on the 
6= 45° meridian of the body of the combination in figures 
35 (o) and 36 (f). The effects shown by the figure are con- 
sistent with those shown for the angle-of-attack case and 
for the wing-incidence case on the top meridian of the body. 

Pressure distribution on wing of wing-body combination. — 
A comparison between linear theory and experiment for the 
pressure distribution along several spanwise stations is made 
in figures 35 (g) to 35 (n). The experimental data (figs. 
35 (k) and 35 ( l )) show that, in general, fiP/i w for the i w = 
— 1.9° data is constant and nearly equal to —2 in front of 
the Mach cone. Behind the Mach cone the theory generally 
predicts values about 5 percent above the experimental data 
for i w =— 1.9°. The higher-order effects due to i w cause 


larger differences between linear theory and experiment for 
i w =—5.7°. The figures show that these differences are 
well predicted by shock-expansion theory. The effects due 
to the influence of the Mach waves are the same as those 
discussed for the angle-of-attack case. There is no effect 
of Mach number evident on the wing of the wing-body com- 
bination other than that predicted by linear theory. 

Span load distribution. — A comparison between the 
theoretical and experimental results for span load distri- 
bution on the wing and body of the combination is mado in 
figure 36 for i w = — 1.9°. The decrease in the span loading 
due to the wing tip was calculated by the method of Buse- 
mann (ref. 20). In part (a) of figure 36, interference from 
both the body and the wing tip is felt between points 1 and 
2, but in part (b) no interference is felt between points 1 and 
2, and the span loading is that of a two-dimensional wing 
alone. 

Figure 36 shows that, in general, the experiment is 5 per- 
cent lower than the linear-theory prediction. Since all 
pressure measurements for the wing-incidence case were 
made for negative values of i w , the experimental values 
used in this figure were obtained by doubling the values 
of fiPjiw obtained for i w = — 1.9° rather than by considering 
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x/pa 

fi) At =1-48, y/a= 1.92 
(j) At=2.00, y/o= 1.92 

Figure 35. — Continued. 

two surfaces as for the angle-of-attack case. Since this 
increases the nonlinear effects of i w rather than m in imizing 
them, only the iw=— 1-9° data (for which the nonlinear 
effects are small) were plotted. However, the present 
method is applicable to the prediction of the net span 
loading for larger values of iw because the nonlinear effects 
on the upper and lower surfaces tend to cancel each other, 
as shown for the angle-of-attack case. 

EFFECT OF REYNOLDS NUMBER 

The primary effect of Reynolds number in this investiga- 
tion was on the body. Reynolds number was found to have 
no significant effect on the pressure distribution on the wing 
of the combination for the range investigated. Figure 37 
shows the boundary -layer condition, as observed in schlieren 
pictures, on top of the body at the point of intersection with 
the Mach wave from the leading edge of the wing-body 
juncture for R— 0.6 and 1.5X10 6 . The transition and 
separation regions shown in figure 37 indicate approximately 
the ranges of a B and iw in which the boundary layer changes 
from laminar to turbulent or separated flow at the Mach 
wave from the wing-body juncture. In laminar and turbu- 
lent regions the flow remains laminar or turbulent across the 



x/po 

(k) M=1.48, y/a= 2.58 

(l) Af=2.00, y/a=2.58 

Figure 36. — Continued. 

Mach wave. Some of the Reynolds number effect shown 
by figure 37 may be due to changes in the turbulence level of 
the wind tunnel. 

It is to be expected that data obtained for several angle 
combinations within any one of the regions shown in figure 37 
would show no significant differences due to viscous effects, 
but that these data would differ from data in other regions. 
For example, for M=1.48 and R= 0.6X10° the data for 
a B = — 2° with 'i w =0° should differ from the data for 
a B = — 6° with - 1 ^= 0 ° because transition occurs at the shock 
wave for the latter case but not for the former. That there 
is a difference is shown in figure 38 where the pressure dis- 
tributions on top of the body for these two conditions aro 
compared. In front of the shock wave the flow is laminar 
for both angles of attack so that there is no difference in the 
two sets of data. However, for a B =— 6° transition occurs 
at the shock wave and the pressure rises as predicted, while 
for a B = — 2° laminar flow persists behind the point at which 
the shock is expected and the pressure rise occurs much later 
than predicted. In fact, the pressure rise does not occur 
until the transition point shown in the figure is reached, and 
then it tends to overshoot. This phenomenon of the delayed 
pressure rise wus observed to occur whenever laminar flow 







QUA8I-CYLINDRICAL THEORY OF WING-BODY 



(m) Af>=1.48, t//a=3.92 
(u) M=2.00, j//a=3.92 

Figure 36. — Concluded. 


persisted beyond the point at which a shock wave from the 
wing was predicted to exist. When the disturbance from 
the wing is an expansion wave, the pressure-coefficient curves 
rise approximately as predicted, regardless of the type of 
boundary layer. The conditions for which this delayed 
pressure rise was observed to occur are shown by the dotted 
areas in figure 37. Two other examples of this phenomenon 
may be seen near points 1 of figures 33 (d) and 35 (d) for 
2°, i w =0° and a B — 0°, ij 7 = — 1.9°, respectively. 

In figure 39, the pressure distributions on top of the body 
are compared for three Reynolds numbers. It is shown that 
data for the two highest Reynolds numbers, 12=1.2 and 
1.6X10°, agree well, while the data for the lowest Reynolds 
number differ from those for the higher Reynolds numbers. 

COMPARISON WITH THEORY FROM OTHER SOURCES 

The three theories for which numerical results are available 
are compared in figure 40. The theory due to Ferrari was 
obtained by cross-plotting from a figure in reference 22 so 
that the curve shown is only approximate. The theoretical 
curve due to Morikawa is obtained from tabulated results 
given in reference 4. The experimental-data region was 
determined by the extreme values obtained for a B — ±2° for 
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0 12 3 4 5 6 

y/a 


(a) Af=1.48 

(b) Af=2.00 

Figure 36. — Span load distributions due to wing incidence; 

1.6X10®. 

Mach numbers 1.48 and 2.00. From this figure it appears 
that either the theory of Morikawa or the present theory can 
he used to predict the pressure distribution in the juncture 
of a wing-body combination. Ferrari’s theory predicts 
values that are somewhat low at the leading edge of the wing, 
but it appears that if numerical results were available beyond 
xjf3a=0.7, they would lie within the experimental range. 
For a more complete comparison of the theories of Ferrari 
and Nielsen, see references 9 and 23. Except for the present 
theory, no numerical results for the pressure distribution on 
the body were available for comparison. 

CONCLUSIONS 

A theory of wing-body interference for supersonic speeds 
has been developed. The theory was applied to the calcu- 
lation of the separate effects of body angle of attack and w ing 
incidence on the pressure distributions acting on a rectan- 
gular wing and body combination. On the basis of compar- 
ison between the theoretical predictions and experimental 




Angle of attack, ag, deg 
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Laminar 
Turbulent 
Delayed pressure rise region 


Transition 

Separation 





(c) (d) 


Wing-incidence angle, i w , deg 


(a) jtf=1.48, 5=0.6X10* (b) Af=1.48, 5=1.5X10’ 

(c) Af=2.00, 5= 0.6X10" (d) M=2.00, 5=1.5X10’ 

Figure 37. — Boundary-layer condition on top meridian of body at 
point of intersection with Mach helix from wing. 
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Figure 38. — Effect of transition position on pressure distribution on 
top meridian of body; M=1.48, 5=0.6X10’, and iV=0°. 
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Figure 39. — Effect of Reynolds number on pressure distribution on 
top meridian of body; M=1.48, a B ~— 6°, and tV=0°. 



Figure 40. — Comparison among several theoretical calculations of 
pressures in the wing-body juncture of the combination; anglo-of- 
attack case. 


measurements, the following conclusions are drawn : 

1. The present theory predicts the pressure distributions 
due to wing incidence about 5 percent high for angles up to 
2°. However, the pressure distribution due to angle of 
attack is predicted about 5 percent low for M= 1 .48 and about 
10 percent low at M=2.00 for angles between +2° and —2°. 

2. Nonlinear effects due to angle of attack and wing- 
incidence angle are large. On the wing the difference from 
linear theory due to nonlinear effects of angle can he pre- 
dicted by shock-expansion theory, except near the wing-body 
juncture for the angle-of-attack case. 

3. Span loading was shown to be predicted within ±10 
percent for both the body and the -wing. The predicted span 
loadings are high for the wing-incidence caso and low for the 
angle-of-attack case. 

4. For the angle-of-attack case, the pressure coefficients 
on the wing are experimentally about 5 percent higher for 
Jlf=2.00 than for M=1.48, when reduced to a form that is 
theoretically independent of Mach number. Otherwise 
Mach number has no important effect. 

5. Viscous effects are important only on the body where 
the shock wave from the wing causes large boundary-layer 
and shock-wave interactions for some angle conditions. 
Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Jan. 4, 1954 







APPENDIX A 

DECOMPOSITION OP BOUND ABY CONDITIONS OF WING-BODY COMBINATION 


A detailed analysis of the boundary conditions for a wing- 
body combination is now carried out for the following 
conditions : 

1. The wing is a flat plate in the 2=0 plane. 

2. The body is an infinite cylinder, the r=l cylinder. 

3. The leading edges of the wing are supersonic. 

4. The Mach number is -^2. 

Consider a wing-body combination corresponding to figure 
4 (a) and shown in greater detail in figure 41 (a). The 
potentials for the flow must satisfy several conditions: 

1. It must be a solution of the wave equation. 

2. It must produce no flow normal to solid boundaries. 

3. It must produce no upstream-moving disturbances. 8 
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and on the body 

^=0; r-1 (A2) 

The first step in the decomposition is to break <p„ into a 
potential due to flow along the x axis and one along the z 
axis in accordance with the superposition principle 
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(b) Parallel flow. 

(c) Crossflow 

Fiqube 41 . — Decomposition of boundary conditions. 


If ip a is the potential for the complete flow about the wing- 
body combination, then the boundary condition on the 
wing is 


^=-i r V; z=0 


(Al) 


The flow conditions at infinity and the prescribed normal 
velocities at the combination surface also obey the super- 
position principle. 

The next steps in the decomposition are to resolve <p b 
and <p c into potentials that can easily be computed. The 
decomposition of ip b into wing-alone and body-alone problems 
is illustrated in figure 41 (b). 


Potential: 
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The potential <p t 

due to the 

wing 

alone at incidence 


produces a velocity field normal to the r=l surface to bo 
occupied by the body. The normal velocity field is decom- 
posed into a Fourier series. Since the wing leading edges 
are supersonic, we can consider the flow above the 2=0 
plane alone. To preserve the wing-alone boundary condition 
when tp d is added, we must confine ourselves to cosine terms, 
and, because of a vertical plane of symmetry, we must retain 
only cosine terms of even multiples of 6 in the Fourier series. 
To counteract the distortion of the r=l surface due to the 
wing alone, a body with opposite distortion is added in the 
form of <p d . 


* Since the surfaces on which tho boundary conditions are given are parallel to the * axis. 
It Is necessary to havo this condition In order to differentiate upstream from downstream. 
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The decomposition of <p e into three components is convenient; a component <p t associated with cylindrical crossflow, a 
component due to a distorted body alone, and a component <pn due to a twisted wing alone. 


Potential: 

Vc = 

Vi 

+ 

Vi 

+ Vn 


0 = 

0 

+ 

0 

+ 0 

(ti)j 

a 3 V= 

* B V 

+ 

0 

+ 0 

(l£L* on ™ g 

0 = 

“* y ( 1+ 1?) + 

0 

- r«,(i+£) 

031 °y lill<ier: 

0 = 

0 

03 CO 

+~!>~) fi.(x) cos 2710— T) F.(x) cos 2 n9 

n— o n—o 


The crossflow associated with causes an upwash distribu- 
tion in the z=0 plane which requires an equal and opposite 
twisted wing to counteract it. Again the r=l surface is 
distorted by the wing-alone flow field, and a cylinder with 
opposite distortion is introduced in the form of <p„ to counter- 
act the distortion. 


Two convenient cases in the wing-body interference 
problem are differentiated; the wing-incidence case in which 
a.B=0 but and the angle-of-attack case in which 

ijr = 0 but oijijAO. The wing-incidence case is represented 
by <pn and the angle-of-attack case by <p c . 


APPENDIX B 


ASYMPTOTIC SPAN LOADING FOR FIRST HARMONIC; WING-INCIDENCE CASE 


The span loading for the first harmonic as given by 
equation (57) is 


These integrals are given in terms of Anger and Weber func- 
tions as given on page 310 of reference 24. 




(Bl) 


Let us take a—1, p—1, and P n =—P L =P„] then the span 
loading is 



J sin 0 cos (is sin 9) dd—trEi (is) 
J sin 0 sin (is sin 0) dO=irJi(is) 


(Weber function! 


(Anger function) 


The value of F 0 (s) is then 


The method of calculating the asymptotic formula for the 
span loading is first to expand the Laplace transform in a 
series about the origin and then to take the inverse transform 
term by term. 

From equation (17) 

«> 

It is now necessary to find the series about the origin for the 
two parts of the transform of equation (B3) . 

From equation (47) there is obtained 
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o 

For small values of s the Anger and Weber functions have 
the expansion 


Ei (m)«- [1+0 (is) 3 ] 

TV 
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so that 

|+0(s 3 )] (BO) 

The ratio of Bessel functions given in equation (B3) has 
an expansion around the origin which is a doubly infinite 
series of products of powers of s and log s. 
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Taking the inverse transform, with the help of reference 25, 
page 282, we obtain the desired result. 


QUASI-CYLINDRICAL theory of wing-body interference at supersonic speeds 
From equations (B9) and (BIO) we get 
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-0(s*)^J [s(y-flogr/2)+ 
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For any Mach number and body radius 

-V 2 
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2xj Pa J 1 (x/j9a) 83 pa 
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APPENDIX C 


DETERMINATION OF WING-ALONE POTENTIAL FOR ANGLE-OF-ATTACK CASE 


The first stop in calculating the potential for the wing 
alone will be to sot up a mathematical model. Since the 
exposed wing of the combination operates in the body upwash 
field which effectively twists the wing, the wing-alone model 9 
is considered to be twisted in the manner predicted by Beskin 
upwash theory for y>a 


aw— a B (l -{~pj—a B ^1 (Cl) 

The concealed wing may bo extended through the body 
region in any manner but, since equation (Cl) gives a w =2 a B 
at both wing-body junctures, it is taken as a flat plate at 
angle of attack 2 a B (see fig. 42). 



Figure 42. — Shape of wing alone with effective twist produced by 
body upwash field; angle-of-attack case. 

Tho twisting of the wing is accomplished by superimposing 
a series of flat-plate wings upon a basic flat-plate wing at 
a w = 2a B (see fig. 43). Each of the superimposed wings is 
at an incremental angle of attack, and each successive wing 
terminates at a value of y greater than the previous one. 
As those incremental values of a become infinitesimally small, 
the resulting potential approaches that of a wing with the 
twist defined by equation (Cl). 

For the purposes of determining the wing-alone potential, 
tho wing is considered to be composed of the three parts 
shown in figure 42: the right exposed half-wing, the left 
exposed half-wing, and the wing section inside the body. 
The perturbation velocity potential is determined for each 


of these wing sections and the results added together to 
obtain the potential for the entire wing alone. Thus, 


<Pw — Vw t + <Pw R 


(02) 


Since the wing may bo considered to be composed of an 
infinite number of flat, rectangular wings, the expression 


<p(.aw,y)- 


■*?(- 


-x cos 


-i ~V 




-y cosh -1 


Vy*+2* 


s cos 


—xy 




=») 


(C3) 


from reference 6 for the velocity potential of a fiat, rectan- 
gular wing will be used as the basic relation for the calcula- 
tions. Equation (C3) gives the velocity potential at any 
point (a* ,y,z) due to a flat, rectangular wing at angle of 
attack <x w , terminating at y= 0, and extending to °° along 
the positive y axis. 10 

Since the twisted wing was shown to be equivalent to a 
basic flat-plate wing at angle of attack 2 a B plus an infinite 
number of modifying flab-plate wings (see fig. 43), the po- 
tential of the right half-wing may be written as 


N 


<Pw B =<p(2a B ,y—l)+J2 v(^a t ,y—T] t ) 


(04) 







ci(y- 2 ag 


Basic wing 


a H' = 2a s + Aa|+Aa2 — 7 
/ 

/ 

iVi^ / ' 



I L -a w *Za B 


2 dg + Ad| 


toJ 

Twisted wing 


Figure 43. — Formation of twisted wing by superposition of infinite 
number of flat plates. 


* Both a and p ore taken as unity. 


The top and bottom wlng^suriaces are still considered independent. 
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Tlie first term on the right in equation (C4) is the potential 
due to the basic exposed half-wing which terminates at the 
wing-body juncture, y= 1, and is at angle of attack 2 a B . 
The second term is the potential of the N modifying wings 
each at angle of attack A a { and terminating at y=y t where 
Since equation (C3) is homogeneous with re- 
spect to a, equation (C4) may be written 

N 

< pw R = , f'& a B,y— i)+S ,y—yi)ha t (C5) 

From equation (Cl) 

da.w= — dy (C6) 

Therefore, 

<Pw R —<p(^aB,y—l)~^a B ^ U <p(l,y—y) ~ 

where the limits of integration are determined by the range 
of y on the wing included in the fore Mach cone originating 
from the point for which <p WR is being determined. From 
figure 44 it is apparent from the equation of the fore Mach 


z 



Figure 44. — Upper value of y included in the foro Mach cono omonat 
ing from Pi) /9=1. 

cone emanating from the arbitrary point P, for /3=1 that 


y a =yi+-yjx 1 i — 2i ! (08' 

The lower limit of integration is at the wing-body juncture 
y= 1. From equations (C3), (C7), and (C8), 


(C7) 


xi 2 =zi s +(v«—yi)* 

Therefore, the upper limit of integration is 


Vwj , — P-P . 


f2Va 


,r «n 
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-z cos 
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•z cos 
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V(y— vf+ zVP— 2* J 
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( 00 - 


Carrying out the integration and combining terms gives 


plane, the potential for the other half-wing is simply obtainec 
by replacing y by —y in equation (CIO). 

<Pwjrr=<Pw R (x,-y,z) (C12; 


(~Sj-3~y) cosh -1 -j= X ■ ■ = ■ 
V VCy-D’+z* 


Combining equations (CIO), (Cll), and (Cl2) gives the 
potential due to the entire wing alone. 


% cosU-i 

(y’+Ph/P-^-z 2 


( z+ y-+z) 


cos 


— x{y— 1) 


V(y— i) j +pVp 




(CIO) 


- BP - -V s {* ( cos '' v& +co,r 1 ift)+ 

(tf+?+y- 2 )cosh- 


V(y-1) 3 +P 


Equation (CIO) gives the potential due to the exposed right 
half-wing. To this must be added the potentials due to the 
other two wing sections. The potential due to the section of 
the wing in the body region is simply the difference between 
the potentials of two flat wings at a w = 2a B - One of these 
wings terminates at the wing-body juncture at y= — 1, 
figure 42, and extends (through the body) indefinitely in the 
positive direction. The other wing terminates at the other 
wing-body juncture, y=+l, and also extends indefinitely 
in the positive y direction. The difference between the 
potentials of these two wings is the potential of the wing 
section in the body 


‘Pw B =<p(2 a B ,y+l)— <f>(2a B ,y— 1) (Cll) 

The expression for <p {a w , y) is given by equation (C3). 
Since the model is symmetrical about the vertical y=0 
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cosh 


(018 


Investigation of ^tr as given by equation (013) rovealf 
that there are three regions on the body in which the roa. 
part of this expression assumes different forms. A fourth 
region, region IV, is entirely on the wing and is, therefore 
not necessary for determining the normal velocity distribu- 
tion on the body. These regions are determined by three 
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FiQtmB 46. — Intersootion of aharaoteristio Mach cones with *=0 plane 
showing corresponding regions. 


characteristic Mach cones. One of the Mach cones origi- 
nates on the body axis at x=0, and the other two originate 
at the leading edges of the two wing-body junctures (see 
fig. 46) . The expressions for the real part of <p w in the three 
regions on the body are: 

Region I : 

w ~ J +coa " JS§] - 

(■ - j +y— 2) cosh -1 ... X ==+ 

Kif+z 2 J J Vfy-i^+z 2 
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Chart 1 . — Values of W Jn ( x , r ) functions. 
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bSbbibbibbibbbbibbbbbbbbSbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbI 

■■iiiiiiiiiiiaBiiaiii»iiiiBiBiBaaaai»iiBiiiaiiBBia»iiii»iaiaaaiiiiiiaiaBiii»iiaiaa»i»iial 

■■bbibmibbbbbbbibbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbmbbbbbbbbbbbbbbbbbbbbbbbbbbbbbI 

bSbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbibbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbI 

■ BiaiBfllBBUBBBIBBBBiaBBBIBBaBBBBBBBBBBBBflBBIBBBaBBflBBBBBBBBaBIflflflBBBflBBBUBBBBBBBBBBBBBBBBflaBBl 

b5bbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbb| 


I ia ■■■■■■■■■■■■■■ ■■■■BBBBimmimmimmimii ■■■■■■■■■■■ 

IbbbbSbbbbbbbbbbbbbbbbbbbmbbbbbbbbbbbbbbbbbbbbbbibbbbbbbbbbbbbbbbbbbbbbbbbibbbbbbbbbbibbbbbbibshI 

l»BBIBBflBBBflBflBflflBBflBBflBBBBflBBBflBBBBIB 


3.6 3.8 


x < 4 
ued. 




W z [x,r) 
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(c) TFj(x,r) 0 < x < 2 
Chart 1 . — Continued. 



Wo(x.r) 



(d) W 3 (x,r) 2 < x < 4 
Chart 1 . — Continued. 



W 4 (*, r ) 
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(e) Wi(x,r) 0 < x < 2 
Chakt 1. — Continued. 


^4 U,r) 
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WAx.r) 
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TABLE n— VELOCITY AMPLITUDE FUNCTIONS 



413072—57 85 
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TABLE ni— ORIFICE LOCATIONS ON WING AND BODY OF WING-BODY COMBINATION 


[Dimensions in inches measured from wing leading edge] 


Vla= 1.020 0 

yfa= 1.253 

y/a= 1.916 

y/a= 2.583 

j//a=3.916 

0=90° 

0=46° 

0. 400 

0. 400 


0. 400 


-0. 281 

0. 408 

.775 

.775 

0.775 

.775 


. 219 


L 150 

L 150 

L 150 

L 150 

L 150 

. 719 

1. 408 

L 525 

L 525 

L 525 

L 525 

L 525 

1. 219 

1. 008 

L 900 

L 900 

L 900 

L 900 

L 900 

L 719 

2. 408 

2. 275 

2.275 

2. 275 

2. 275 

2. 275 

2. 219 

2.008 

2. 650 

2. 650 

2. 650 

2. 650 

2. 650 

2.719 
3. 219 

3.719 



















TABLE IV.— PRESSURE COEFFICIENTS 

(a) AT-U8, B-iJSXio> 
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TABLE IV. — PRESS URE COEFFICIENTS Continued 

(b) «■=.!.«, Ji^uxioo 



1352 REPORT 1262 — NATIONAL ADVI80EY COMMITTEE FOR AERONAUTICS 










































































TABLE IV.— PRESSURE 


(o) AT- 1.48, 


8t»ll0O 

\ hi 
*\ 



-0* 




— 1* 




— 2° 

-5,7* 

-3^“ 

-I.9* 

0° 

-6 .r 

— 38* 

-1.9* 


69 

-3^* 

—IIP 


*/»*< 0.170- — 

0. 138 

0.0© 

mm 

1431 



EPfi! 


a 9i« 

a. 390 

RjS 




,158 

.cos 

.406 

.324 

.255 

.578 

■EE 

PBPT7 

. 173 

ISBImS 

.m 

.14? 



• SSI 

^■rTTT 

.374 

.204 

.£24 

.344 

.»* 

■CmTT 

.148 

^K'T? 

,177 

. US 



.w$ 

.450 

.174 

.308 

.221 


.977 

.1©$ 

.144 

,133 

.171 

.114 



.CR 

.444 

.374 

.112 

. 549 

.346 

.291 


.140 

.260 

■ A] 

.U1 



■ TO 

• 4W 

■Ik 


.251 

.342 

• »7 


.172 

. 261 

.20! 

, 124 



.881 

.383? 

.313 

.276 

.227 

■yj 

. 950 


.103 

.131 

.186 

.133 


jr/*-0J08 

.135 


.401 

.388 

.316 

.446 

.358 

.176 

. .210 

m 

.236 

.15© 



r* 

.510 

.12$ 

. 33$ 

.■m 

,3©6 

.311 

.141 

.183 

.273 

.216 

.14© 



.3N 

.478 

.399 

.314 

.243 

.367 

.9* 

.316 


.247 

.101 

.133 



.see 

.461 

. 378 

.316 

. 237 

.3M 

.237 

.206 

.IK 

.247 

.17! 

.111 




.439 

.371 

.312 

.240 

■EE 

.980 

.217 

.145 

.263 

.186 

.118 



. 758 


.JO 

.20© 

. 246 

.837 

. .270 


^■vTT 

. 25C 

.105 

.123 



.m 

.*7 

.339 

.271 

. 213 

,312 

.261 


.3t& 

.232 

.177 

.116 




.956 

mm 

ms 

n 

,281 

.427 

.314 

.m 

.101 


.933 

.167 




lE 

.407 

.348 

.5B9 

: {^WTir 

.821 

.214 

s^BtvTT 

.988 

.914 

. 144 



.ore 

lu 

.171 


. 260 

.M3 


.21$ 

.167 

. 976 

.213 

. 117 



MMMmlr?; 


.144 

.*9 

.340 

.527 

gB£77i 

.224 

.188 

.2X8 


. 106 



.75© 

.195 

.321 

■ 1541 

■E£ 

.903 

. 241 

. 185 

.148 


. 165 

MBS 



.884 

.179 

.307 

■n 

.168- 

.283! 

* 935 

.167 

.194 

. 183j 

.143 

Mm 


n/4-auu 

.131 

.641 

. 651 

.417 

.2971 

. 516 

.384j 

.261 

.176 

.344 

.245 

.lari 



.158 

.643 

.484 

K 

■fTTi 

.464 

.371 

. 274 

.197 

.343 

. 265 

. 1M 



.383 

.47© 

■TP7 

.171 


. 408 

.3*9 

.273 

.1E£ 

■Kvil 


. i as 




.430 

.361 

.51© 

.172 

. 343 


.958 

.106 

.987 

liBFTT 

. 16S 



.631 

.180 

■EE 

.355 

.213 

.281 

.m 

^B-i v 

.16© 

.118 

.1M 

, 135 



.75© 

.MB 

. 285 

.126 

.178 


PHt/j!-. 

.163 

.117 

,177 

.147 

B 



. txn 

. cg« 

. 387 

.221 

. x\jw 

. *Oij 


|A4» 

. 110 

. 168 

. 1A» 



t/lmOMa 

.333 

.433 

H 

.123 

. 157 

. 558 

.311 

■RT 

.163 

. 280 

.21© 

BEE 



.506 

.373 

.316 

. 272 

.124 

^■FjV 

Mra 

■STf 

,1M 

.939 

,1W 

■fETi 



.633 

.343 

.272 

.236 

.188 

.257 

.117 

i^BKTo 

.131 

,19# 


.1® 



.70 

.334 

.216 

.197 

.1® 

.199 

.188 

{■fcrt 

.110 

,168 

• 13ff 

.009 



.684 

.31© 

.241 

.181 

» 141 


.187 

■iLii 

.098 

,151 

.11© 



9-80* 

-.004 

■fiTff 

.018 

■ ' 

mbp 


.011 

Hliiv 

ms 

.01© 

.Oil 

.006 



.073 

LyQiTJ 


.061 



. 066 


■E3 

.053 





■ HO 

.375 

■BE 

.067 



.074 

^BCTTiri 

.DM 

.050 

.056 

.048 



.410 

.38© 


.182 

,0S8 

.991 

, 172 

.074 

.006 

.144 

. W£ 

.049 



.571 

.43© 

.383 

.287 

■Lu 

.335 

.975 

.186 

.001 

.24 1 

.186 

, 040 




.46© 

. 367 

.338 

.380 


.3215 

.237 

.107 

.117 

.216 

.133 



. 800 


.381 


.377 

.m 

.313 

.267 

.131 

.102 

.240 

.154 



L 073 

.343 


.«e 


.995 

.262 

.224 

.396 

.246 

.219 

.173 



1.140 

.172 

.178 

. is© 

.178 

. 150 

.154 

.164 

.146 


. 126 

.US 


l“«* 

.166 

.133 

■7 

.063 

.067 


.085 

■P5 


■53 

■ 5 ? 

■PFp 




. 421 

.»7 

.177 

. 07© 

.270 

.106 

.078 

.049 

.100 

.086 

.0® 



.48© 

.471 

.408 


.115 

■EE 

.271 

.190 

. 158 

^■£v_P 

.179 

.110 



.656 

.451 

.405 

.133 

. 148 

.380 

. 315 

,ai 

.118 

.278 

.181 

. 100 



.83 

. 4U 

.374 

.334 

. 173 


. 999 


.184 

.973 

.2*7 

.141 



.990 

.344 

.115 

.174 


,2M 

.260 

, ©IQ 

.109 

.930 

.199 

.141 
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